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ABSTRACT 



We present an optical/near IR selected catalog of 79 clusters distributed over 
an area of 5.1 square degrees. The catalog was constructed from images obtained 
with the 4-Shooter CCD mosaic camera on the Hale 5m telescope operating in 
"scan" mode. The survey, hereafter known as the Palomar Distant Cluster Survey 
(PDCS), was conducted in two broad band filters that closely resemble V and /. 
The 4(7 limiting magnitudes for our 300 s exposures are ~ 23.8 (V) and 22.5 (/). 
A matched filter algorithm was developed and employed to identify the cluster 
candidates by using positional and photometric data simultaneously. The clusters 
cover the range 0.2 ^ z ^ 1.2, based on the redshift estimates derived in the 
cluster detection procedure. An accurate selection function is generated from 
extensive simulations. We find that the cumulative surface density of clusters 
with richness class R > 1 is about a factor of 5 higher than the extrapolated 
density of R > 1 Abell clusters. The PDCS results are consistent with a constant 
comoving density of clusters to z ^ 0.6, albeit at the above high density level. 
Constraints on cluster abundances at 2; > 0.6 will be possible with the acquisition 
of spectroscopic redshifts for a large subset of these cluster candidates. We also 
present a supplemental catalog of 28 clusters that do not satisfy all our selection 
criteria but which include some of the most distant systems detected in the 
survey. Finding charts for all 107 cluster candidates are provided. (Astro-ph 
readers: due to their large size, the finding charts are not included in this version 
of the paper). 





Subject headings: galaxies: clustering; cosmology: observations 
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1. INTRODUCTION 



The evolution of the space density of clusters, of their internal kinematics and spatial 
structure, of their large scale clustering properties, and of their host galaxy populations are 
fundamental predictions of all viable galaxy and cluster formation theories. The accuracy of 
the observational constraints on these parameters depends on the availability of statistically 
complete cluster catalogs over a large redshift range (0 < 2; ^ 1 and beyond). At redshifts 
less than 0.2, the available catalogs (Abell 1958; Zwicky et al. 1968; Shectman 1985; Abell, 
Olowin, & Corwin 1989, hereafter AGO; Lumsden et al. 1992 (Edinburgh-Durham); Dalton 
et al. 1994 (APM)) have all been based on the large photographic sky surveys done with the 
1.2m Schmidt at Palomar Mountain and at the UK Schmidt at Siding Spring. These surveys 
have only recently become available in digital format and, consequently, only the Shectman, 
APM, and Edinburgh-Durham catalogs have been generated using automated and objective 
cluster detection algorithms. At higher redshifts, the catalogs are based on deep prime 
focus optical/near IR photographic plates taken on 4-5m class telescopes (Gunn, Hoessel, & 
Oke 1986, hereafter GHO; Gouch et al. 1991). The low quantum efficiency of photographic 
material limits the completeness of these surveys beyond z ~ 0.5. These catalogs have 
yielded important information on galaxy evolution {e.g., Butcher & Oemler 1984; Gunn & 
Dressier 1988; Rakos & Schombert 1995) and on the clustering properties and space densities 
of low redshift clusters as a function of richness {e.g., Bahcall & Soneira 1983; Dalton et al. 
1992; Nichol et al. 1992; Postman, Huchra, & Geller 1992). The Abell/AGO and Zwicky et 
al. cluster catalogs and the previously published optically selected cluster catalogs at higher 
{z ^ 0.3) redshifts were constructed via visual inspection, limiting their utility. Eor a cluster 
catalog to be statistically useful, it must be homogeneous and its selection biases must be 
accurately quantified. Eurthermore, in order to compare the predictions of models to the 
observations reliably, the same cluster selection process that is used on the observations 
should be applicable to large scale simulations. 

The 4-Shooter GGD mosaic camera on the Hale 5m (installed in 1983) was designed, in 
part, to conduct digital, two color surveys for distant {z ^ 1) clusters over a moderate area 
of sky. We present here the cluster catalog resulting from the first such survey, hereafter 
referred to as the Palomar distant cluster survey (PDGS). We have developed a matched filter 
cluster selection algorithm which uses both positional and photometric data simultaneously 
to search for overdensities in the galaxy distribution. The algorithm works well over a large 
dynamic range in redshift and is optimal in the regime where one is looking for small signals 
in large backgrounds, as is the case for distant clusters. We have applied this algorithm to 
the uniform galaxy catalogs derived from our GGD data. 

Section 2 of this paper describes the observations and their calibration. The construction 
of galaxy catalogs from the GGD data is discussed in section 3. Our cluster selection 
algorithm, the cluster selection criteria, the cluster catalog, the catalog completeness and 
the false positive rate are discussed in full detail in section 4. Preliminary constraints on the 
space density of optically selected clusters are discussed in section 5. The properties of the 
clusters in this catalog will be discussed in subsequent papers. 
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2. OBSERVATIONS AND DATA CALIBRATION 

2.1. Observational Technique 

The PDCS consists of five separate fields. The shape of each field is roughly a square 
that is f degree on a side. The fields are all located within or near the original high-latitude 
Gunn & Oke (f975) survey areas. In addition, field centers were chosen to minimize the 
number of bright stars in the survey regions and to avoid very low-redshift [z < 0.1) clusters 
and regions with high reddening. There is one R = f Abell cluster, A 1739 {zest ~ 0.17), 
located in our survey area. The coordinates of the field centers are given in Table I; the 
five regions will be denoted in this paper as the 0^, 2^, 9^, 13^, and 16^ fields. The 2^ CCD 
survey field is not included in the original Gunn & Oke survey; this field is used in lieu of 
the original 3^ field which suffers from high galactic extinction. 

The PDCS was performed with the 4-Shooter CCD camera (Gunn et al. 1987) mounted 
at the cassegrain focus of the Hale telescope. The focal plane of the telescope is coincident 
with a refiective four-faceted pyramid in the 4-Shooter. The pyramid directs the incoming 
light into four reimaging cameras. The 4-Shooter detectors are four 800 X 800 Texas 
Instruments CCDs. The cameras produce an image scale of approximately 0.335" pixel"^; 
the instrument's entire field of view is 8.9' on a side. A diagram of the 4-Shooter focal plane 
is given in Figure I. 

The CCDs are positioned so that there should be no gaps in sky coverage along the 
pyramid facets (there are typically about 20 rows/columns of pixels that do not see the sky 
due to pyramid shadowing), but the few arcseconds of sky located where the CCDs join are 
not usable because of imperfections in the pyramid surface along the facets. The effective 
size of the field for this project is 2 X 780 pixels, or a width of 520". Images of each field 
are obtained in two filters, the F555W and F785LP of HST's Wide Field/Planetary Camera 
(Griffiths 1990). To minimize confusion with the HST photometric systems, we will denote 
the 4-Shooter passbands as (F555W) and /4 (F785LP); see §2.2 for a description of the 
PDCS photometric system. The limiting magnitude of the data is fainter than that of the 
/4 observations (due to the fainter sky and higher detector sensitivity at 5500 A), but the /4 
images are in some respects just as important; ai z ^ 0.6 the 4000 A "break" lies redward 
of the band, so the /4 data provides the wavelength coverage needed to detect clusters 
out to 2; ~ I. 

The data are acquired in Time-Delay-and- Integrate, or TDI, mode. (This is frequently 
called "scan" mode in the literature.) Details of how TDI observations are taken with the 
4-Shooter, as well as the initial data processing of the data, are given in Schneider, Schmidt, 
& Gunn (1994). In TDI mode the sky is moved along the columns of the CCD at a constant 
rate; the CCDs are read so that the charge packet shifts one row in the time that it takes the 
sky to traverse one pixel. TDI mode was chosen for the PDCS for three reasons: I) This is a 
very efficient method to image large sections of sky to moderate depth (there is no overhead 
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associated with CCD readout; one is exposing essentially 100% of the time), 2) Data taken in 
TDI mode can be extremely well flat-fielded (see Schneider, Schmidt, & Gunn 1994 and §2.3 
of this paper). This property is especially important for these survey images, as the Texas 
Instruments CCDs have very poor fringing properties at near-infrared wavelengths. 3) The 
configuration of the 4-Shooter focal plane (Figure 1) is well-suited for two-filter observations 
in TDI mode. 

The PDCS data were obtained by moving the telescope north along lines of constant 
Right Ascension (equinox 1950.0) at a rate of 0''893 s~^; this yields a CCD transit time of 300 
seconds. Because of the aforementioned pyramid shadowing in the 4-Shooter, the effective 
exposure time is 293 seconds. The data are acquired in two colors, nearly simultaneously, 
by making use of the fact that the 4-Shooter can be equipped with one set of filters for the 
leading two cameras (CCDs 1 and 4, see Figure 1) in the scan and another set for the trailing 
ones (CCDs 2 and 3). 

A survey field consists of eight adjacent scans. Fach scan requires approximately 
80 minutes to complete and overlaps the adjacent scan by 30". The total area imaged in 
each field is usually slightly larger than a square degree (see Table 1). 

Imaging of the five fields began in October 1986 and required approximately three years 
to complete. The data were all obtained under photometric conditions. In order to maintain 
accurate star-galaxy classification and uniform completeness limits, only data with FWHM 
seeing better than 1''5 are used in this analysis. The median FWHM seeing in all slow-scan 
data obtained during the course of the observations is 1''40; the FWHM median seeing in 
data accepted for use in our cluster survey is 1''23. 

2.2. The PDCS Photometric System 

The system response as a function of wavelength for the filter-|-4-Shooter-|-telescope 
are shown in Figure 2. The lower curve for each filter is the quantum efficiency through 
1.2 airmasses from Palomar; this is the response curve used to calculate filter properties 
{e.g., effective wavelengths) and transformations to other photometric systems. 

The effective wavelengths and widths of the filters, following the definitions in Schneider, 
Gunn, & Hoessel (1983) and Grifliths (1990), are 5430 A and 1230 A for V4 and 8690 A and 
1400 A for I4. The effective wavelength of the V4 band is 100 A bluer than that of and 
V4 is about 50% wider than V. I4 has nearly the same width as the Kron-Cousins / band 
(Cousins 1976; this band will be denoted as /kc), but the effective wavelength of I4 is about 
500 A redder than that of /kc- The I4 band has some sensitivity at wavelengths longward 
of 10,500 A. 

The zero points of the V4,/4 magnitudes are based on the AB magnitude system of 
Oke & Gunn (1983). Objects whose spectral energy distribution are constant in terms of 
ergs s~^ cm~^ Hz~^ have neutral color in an AB-hased system. The magnitudes of Vega 
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are V4 = +0.03 and I4 = +0.46. The V4,/4 magnitudes for the fundamental standard of the 
Gunn system, the F subdwarf BD +17 4708 (defined as magnitude 9.500 in all bands), are 
9.45 and 9.20, respectively. Magnitudes in the V4,/4 system for the six standard stars that 
provided the photometric calibration for the PDCS are listed in Table 2. (Kent 1985 gives 
coordinates and g^r magnitudes for these stars.) 

Since the V4 and I4 filters are quite broad, it is impossible to define precise 
transformations to other photometric systems when dealing with galaxies with a wide range of 
redshifts and spectral energy distributions. We have calculated approximate transformations 
from V4 and I4 measurements to a number of other photometric systems that are frequently 
used in studies of faint galaxies and/or clusters. In addition to V and /kc, we present the 
transformations to Bj (Agfr = 4390 A, GuUixson et al. 1994), and the g (Agfr = 4910 A), 
(Agfr = 6530 A), and i (Agfj = 7860 A), bands of Thuan & Gunn (1976) and Wade et 
al. (1979). (Also see Kent 1985 for further information on the g and r photometric 
systems.) The rough transformations were derived by synthesizing magnitudes from the 
spectrophotometric atlas of stars of Gunn & Stryker (1983), and restricting the fits to objects 
with {V4 -l4< +3). 



Bj = V4 + 0.29 + 0.676 {V4 - h) - 0.083 {V4 - hf (1) 

g = V4 + 0.04 + 0.295 {V4 - I4) (2) 

V = V4 - 0.02 - 0.056 {V4 - I4) + 0.012 {V4 - hf (3) 

r = V4 + 0.13 - 0.514(1^4-/4) + 0.055(1^4-/4)^ (4) 

I = h + 0.27 + 0.145(^4-/4) (5) 

/kc = h - 0.43 + 0.089 {V4 - h) (6) 

The average rms scatter about the above relations is a few hundredths of a magnitude. The 
correction for Galactic extinction for the V4 and h filters is 

Av, = Av = 3.05 E{B - V) (7) 

Aj, = 0.497 Av, = 1.52 E{B - V) (8) 



where E[B — V) values are derived from the procedure presented in Burstein & Heiles (1982). 
Reddening corrections for the five survey fields are given in Table 1. 

Figures 3a-c show the V4 and h k-corrections and expected V4— /4 colors as functions of 
redshift, respectively, for five morphological types: F, Sab, Sbc, Scd, and Im. The spectral 
energy distributions are from Coleman et al. (1980) and no spectral evolution is assumed. 



2.3. TDI Image Data Processing 



The raw data from a 4-Shooter TDI observation form a continuous pixel stream 816 
columns wide and approximately 51,000 rows long (800 rows are generated every 300 seconds 
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for each CCD). To take advantage of the large amount of existing picture-based image 
processing software, the data from each CCD are subdivided into 816 X 800 frames with 
an 80 pixel (26''8) overlap introduced in the scan direction. The overlap allows reliable 
photometric and astrometric measurements of essentially all of the objects that fall near 
the boundaries of the image frames. The size of the overlap was chosen as a compromise 
between processing efficiency (larger overlaps mean more redundant processing) and the 
typical diameters of galaxies in moderate to high redshift clusters 15"). Analysis software 
may not accurately measure the parameters of galaxies with diameters larger than 80 pixels 
that are located in the overlap area. The number of such cases is extremely small (< 10 per 
square degree), and in any case they are not the targets of study of this program. 

In TDI mode, the data can be used directly to construct a one-dimensional "skyflat" 
vector that will produce extremely accurate measurements of the variation in pixel-to-pixel 
sensitivity perpendicular to the scan direction. No adjustment for pixel-to-pixel sensitivity 
variations parallel to the scan direction is required in TDI mode. The median sky values 
in the V4 and I4 bands are 2400 DN pixel"^ (~21.3 mag per arcsec^) and 6000 DN pixel"^ 
(~19.3 mag per arcsec^), respectively; thus there is sufficient signal to generate accurate 
flat-field vectors. 

To correct for the relative zeropoint offsets between the two CCDs assigned to a given 
passband (e.^., CCDs 1 and 4 for the V4 band), a sequence of raw frames from each CCD are 
debiased and then fiattened by multiplying the data by their corresponding skyfiat vector. 
A "balance factor" for each CCD, BFi^ is computed using the expression 



I Ef =i(Sky,, + Sky 
Ef=i Sky, 



BF, = ^^"=^^^-^"'' ^ (9) 



where Skyk,i is the sky level in the k^^ frame from CCD z, Sky^^ is the sky level in the k^^ 
frame from CCD j, where j is the CCD assigned the same filter as CCD z, and N is the total 
number of frames processed for this calibration (usually four frames per CCD). The initial 
skyfiat vector is then multiplied by the balance factor to produce the final, corrected skyfiat 
vector. This method of computing relative offsets distributes the offset correction between the 
two CCDs as opposed to applying the correction all to one CCD. The assumption underlying 
this approach is that sky level differences between corresponding frames from CCDs with a 
given passband are due to different relative CCD zeropoints rather than actual sky variations 
on the scale of a few arcminutes. This quite reasonable assumption is supported both by the 
insensitivity of the derived balance factor values to the specific frames used in the calculation 
and by the excellent agreement between the balance factors derived from the scans and 
those derived from V4 and I4 dome fiat exposures. The typical zeropoint offsets are 0™36 
and 0™12 for the V4 and I4 CCDs, respectively (i.e., CCD 4 is 1.39 times more sensitive than 
CCD 1 in the V4 band, and CCD 3 is 1.12 times more sensitive than CCD 2 in the I4 band). 

The CCD signature is removed from the raw data and, as a final step, the median sky 
level as a function of row (scan direction) is fit to a cubic polynomial. The appropriate linear, 
quadratic, and cubic components are subtracted from each image row. This process removes 
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any medium-term temporal variations in sky brightness and produces frames with extremely 
flat sky levels (sky brightness levels measured in different regions in any given frame agree to 
within 0.002 mag per arcsec^). Most frames are quite flat even before applying the final sky 
correction. The largest temporal variations in sky brightness across a frame seen in the data 
have amplitudes less than 2% of the mean sky level. One consequence of this sky correction 
is that the mean sky level will change in a non-continuous manner from frame to frame. 
Measurements of objects in the 80 pixel overlap region between frames shows that this has 
no infiuence on number counts or object photometry and yields very stable object detection 
performance. 

3. GALAXY CATALOG CONSTRUCTION 

3.1. Object Detection and Classification 

We used a modified version of the Faint Object Classification and Analysis System 
(Jarvis & Tyson f98f, Valdes f982; hereafter FOCAS) to detect, measure, and classify 
objects in the calibrated CCD images. The modifications include new software to perform 
photometric classification of galaxies (discussed later in this section), the addition of 
object catalog parameters related to this photometric classification as well as a new object 
compactness parameter, and a variety of new image statistics, catalog editing, and quality 
assurance programs. The basic FOCAS object detection, point spread function (PSF) 
measurement, deblending, and star-galaxy-noise classification algorithms have not been 
altered. 

As in the image calibration process, the frames are processed individually. The 
object isophotal detection threshold in each frame is set to be ScTgky above sky (typically 
25.7 mag per arcsec^ in the V4 band and 24.8 mag per arcsec^ in the I4 band) and a minimum 
object area requirement of f5 pixels (f .68 arcsec^) is imposed. Because the sky level is lower 
and the system response is higher in the V4 band than in the I4 band, the V4 images are 
deeper. Consequently, about 40% of the objects detected in a typical V4 image are undetected 
by FOCAS in the corresponding I4 image. Nearly all the objects which are only detected 
in the V4 band are faint [V4 ^ 23). About 95% of the objects with V4 < 23, however, are 
detected in both passbands. 

After object detection, merged objects are deblended using the FOCAS splits routine. 
Approximately f5% of all detected objects with V4 < 24: required deblending. Object 
parameters are evaluated (e.^., isophotal, core, and total magnitudes, axial ratio, intensity- 
weighted and unweighted moments, etc.) and the FOCAS resolution object classifier is run 
on all of the objects. The PSF used by the resolution classifier is computed separately for 
each frame to compensate for variations in seeing that can occur during the course of a scan. 
Occasionally there is an insufficient number of suitable PSF stars in a frame for an accurate 
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PSF to be measured; in such cases, the PSF from an adjacent frame (in the scan direction) 
is used for object classification. 

FOCAS classifies objects based on the results of fitting a series of templates derived 
from the observed PSF to each object. Specifically, the templates are parameterized by 

tin) = / X pin/s) + (1 - /) X pin) (10) 

where is the position of pixel z, pivi) is the observed PSF, s is the FOCAS resolution scale 
parameter, and / is the FOCAS fraction parameter. The object classification is determined 
from the values of s and / which yield the best fit to the two-dimensional data. We adopted 
the following classification rules: 



Noise 0.1 < s < 0.7 < / < 1 

Star 0.7 < s < 1.3 < / < 1 

Fuzzy Star 1.3 < s < 10 < / < 0.25 

Galaxy 1.3 < s < 10 0.25 < / < 1 

Diffuse 10 < s < 100 < / < 1 



Practically all "Fuzzy stars" are stellar objects and every diffuse object brighter than V4 = 23 
is a galaxy. The range of the scale parameter used to select stars was chosen to be broad 
enough to account for reasonable PSF variations that occur across a single 4-Shooter CCD 
but narrow enough to minimize galaxy misclassifications. The 4- Shooter re-imaging cameras 
are relatively fast (//1. 84) and the PSF varies in shape from center to edge. In TDI mode, 
PSF variations due to the camera optics are a function of CCD column only. FOCAS 
currently does not allow one to use a variable PSF for object classification on a given CCD 
frame. The scale parameter range used to select stars was, therefore, chosen to be broad 
enough to account for reasonable PSF variations that occur across a single 4-Shooter CCD 
but narrow enough to minimize galaxy misclassifications. We find that there are position 
dependent star classification errors due to the PSF variations but the galaxy classifications 
are not significantly affected. This is primarily because galaxies dramatically outnumber 
stars at the magnitude limit of the present survey. A detailed discussion of the accuracy of 
the object classifications is presented in §3.3. 

Finally, we perform a crude morphological classification of the galaxies based on their 
radial surface brightness profile. For galaxies with azimuthally averaged radial profiles 
extending for 1''34 or more, a comparison is made between the observed profile and a library 
of seeing- convolved reference profiles for exponential disk and de Vaucouleurs (1948) r~i 
models. The reference profiles span a wide range of effective radii and axial ratios. A 
maximum likelihood technique is used to select the best fit model. If the galaxy has an 
axial ratio greater than 2.5, only the disk models are used in the comparisons. A new 
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"total" magnitude is then computed by integrating the profile ol the best fit model. This 
total magnitude, the effective radius of the best fit model, the xl value of the fit, and the 
profile type (disk or de Vaucouleurs) are stored along with the other object parameters. 
The morphological classification algorithm reliably classifies pure disk or pure bulge systems 
down to V4 = 22. At the survey limit, however, only about 50% of simulated r"? galaxies 
are correctly classified. There is a classification bias in the present algorithm that results in 
galaxies with V4 > 22.5 being preferentially classified as exponential disks. 

3.2. Photometric Calibration 

Although all of the data were obtained under photometric conditions, very few 
observations of photometric standards were made on the nights that the scans were 
performed. Consequently, a bootstrap procedure was used for the final photometric 
calibration of the PDCS. 

All eight scans of a given field were placed on a uniform fiux scale by comparing the 
instrumental magnitudes of the objects located in the ~30" overlap region between adjacent 
scans. This procedure is done for both the V4 and I4 data independently. The observation 
dates for some adjacent scans are separated by as much as two years; temporal variations in 
the refiectivity /transmission of the telescope and 4-Shooter optics will alter the photometric 
zero points. If the data are indeed photometric, then a single constant for each scan will 
suffice to place the data on a given instrumental system. If the observations were taken in 
nonphotometric conditions, this offset will vary with time. 

Figure 4 shows the photometric offset as a function of time (for the PDCS observations, 
the declination is linearly related to the time elapsed since the start of the scan) for four pairs 
of adjacent scans. The upper two panels show scans taken on the same night; as expected, 
there is no significant offset between scans. The lower two panels compare data for scans 
taken at times separated by more than a month; while there is a definite change in the fiux 
scale between adjacent scans, there is no evidence for any temporal variations. Objects with 
V4 < 22 were used to measure the offset in the V4 band data [I4 < 21 in the I4 band). The 
observed scatter is consistent with the estimated uncertainties in the FOCAS photometry 
(±0.12 mag for a galaxy with V4 = 22). The handful of points with discrepant offsets 
are from objects which straddle the overlap zones and, hence, their magnitudes cannot be 
accurately compared. These objects are excluded from zeropoint offset computations. The 
typical photometric offsets seen in the data are ib0™09 and ibO™10 in the V4 and I4 passbands, 
respectively. There are some cases where the offsets are as large as ib0™25. There is little 
indication of zeropoint drifts which would occur if the photometric quality of the sky varied 
during the course of a scan. The photometric offsets in each band are normalized to the scan 
with the best seeing; following this adjustment the data in each field are on the same fiux 
scale. 

To place the observations on an absolute fiux scale, we acquired three "fast" scans for 
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each of the five fields. The fast scans were taken in a direction perpendicular to the motion of 
the survey scans. The fast scans have an effective exposure time of 98 seconds, span the entire 
width of the field, and are centered on the northern, central, and southern segments of each 
field. The fast scan data were obtained on three photometric nights (27 September 1987 and 
12 and 13 March 1988) on which we also obtained observations of photometric standards 
(see Table 2 for a list of standards). In addition to providing the absolute photometric 
calibration, the data produced by the fast scans are used to tie the strips together in the 
astrometric calibration (§3.4). 

The accuracy of the photometric calibration of the data can be measured by comparing 
number counts from the five independent survey fields. Figures 5a and 5b show the 
differential galaxy surface density as a function of V4 and I4 magnitude, respectively. Figure 6 
shows the residuals between the mean galaxy surface density for all five fields and the galaxy 
surface density in each individual field as a function of V4 and I4 magnitude, respectively. 
All magnitudes have been corrected for Galactic extinction using the values in Table 1. 
The typical field-to-field number fiuctuations are ±9% over the range 19 < V4 < 23.5 and 
18.5 < I4 < 22. At < 19 {I4 < 18.5), the number of galaxies is too small to yield reliable 
statistics; at > 23.8 {I4 > 22.5) incompleteness sets in. The solid curve shows the observed 
rms fiuctuations averaged over all 5 survey fields. For comparison, the dashed curve shows 
the expected fiuctuations between fields of this area, given the observed angular two point 
correlation function derived from the galaxy survey of Neushaefer, Windhorst, & Dressier 
(1991). The observed fiuctuations are in very good agreement with the clustered distribution 
prediction, indicating that, on average, zeropoint fiuctuations between survey fields must 
be ^ ib0™09. This suggests that the galaxy catalogs are not seriously contaminated by 
misclassified stars, a conclusion confirmed by our checks on classification accuracy described 
in the next section. The systematic excess of galaxies in the 0^ field is most likely due to 
a photometric zeropoint error (there is no particularly rich cluster of bright galaxies in this 
field). A zeropoint error of — 0™20 (roughly a 2a deviation) in both passbands for the 0^ data 
is sufficient to explain the excess at faint magnitudes [V4 > 20) and to make the differences 
at < 20 statistically insignificant. Our cluster identification algorithm is insensitive to 
zeropoint errors of this magnitude. The good agreement in the faint galaxy counts between 
fields demonstrates that on very large scales (^ 600h~^ Mpc) the optical universe does 
appear to be homogeneous. 



3.3. Object Classification Accuracy 



The uniformity of the log N^a/ - magnitude relationship from field to field is certainly 
a necessary but not sufficient condition for confirming good classification accuracy. An 
additional check is to compare object classifications from the V4 and I4 images. We obtain 
identical classifications for 95% of the objects detected in both passbands when the object 
V4 magnitude is 20.5 or brighter. For objects with V4 < 22.5, identical classifications are 
obtained 86% of the time. For the statistically complete survey used for cluster detection 
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[V4 < 23.8), identical classifications are obtained lor 75% ol the objects. 

There are two reasons why object classification accuracy is not perlect. The first is that 
the FOCAS classifier (and, indeed, most any classifier) breaks down lor objects lainter than 
about 1 magnitude above the completeness limit. Near the limit, the signal-to-noise ratio 
is low (S/N ~ 4 at = 23.8) and the mean object area is ^ 5x FWHM^, making reliable 
classification difficult. The second is that the PSF varies in shape from center to edge (c/., 
§3.1). Since we are using a single PSF per frame, classification accuracy should degrade 
somewhat near the east/west edges of every scan (in TDI mode, the PSF is a function of 
column only). We have attempted to minimize the latter effect by allowing a 30% tolerance 
in the fit to the mean PSF for objects classified as stars (c./. §3.1); however, the effect is 
nonetheless still present. Fortunately, the net result is to introduce only minor contamination 
of the galaxy catalogs, as shown in Figures 7 and 8. The star catalogs are more seriously 
contaminated. Statistical corrections are made to compensate for the contamination, as 
discussed below. 

We now demonstrate the nature and amplitude of the misclassifications. Figure 7a and 
7b show the galaxy-to-star ratio (GSR) for the 13^ field as a function of V4 and I4 magnitude, 
respectively (results for the other fields are very similar). The solid points represent the data 
for all objects regardless of CCD column; the open points represent data for objects with 
200 < CCD Column < 600. The breakdown of the classifier at low signal-to-noise levels is 
the source of the turnover in the GSR in Figures 7a and 7b. The probability of a faint star 
being classified as a galaxy is higher than the probability of a faint galaxy being classified as 
a star. This has been verified by running the classifier on simulated CCD images. Galaxies 
outnumber stars by about 7 to 1 at ~ 23 and, consequently, even a small error in galaxy 
classification results in a big change in the faint star counts. For the same reason, the galaxy 
counts are not significantly affected. For V4 ^ 21 the GSR is essentially independent of CCD 
column. At fainter magnitudes, the GSR is systematically lower for objects located within 
the central 400 columns. This is precisely the trend expected from a PSF variation - the 
assumed PSF is rounder than the actual PSF near the east/west edges and hence relatively 
more stars are misclassified as galaxies in those areas. This is shown directly in Figures 8a 
and 8b which depict the number of stars and galaxies as a function CCD column number. 
There is an obvious drop in star counts in the 150 outermost columns. Again, the galaxy 
counts are not as noticeably affected because galaxies vastly outnumber stars at ~ 23 and 
because the broadened PSF is still narrower than the typical galaxy profile. 

Using both the observed data and results from simulations, we have adopted effective 
weighting factors to correct the galaxy counts depending on position and magnitude. At 
V4 = 23.8 and I4 = 22.5 (the limits adopted for the purpose of finding clusters) and near 
the east/west CCD edges the galaxy count correction is typically 15%, only a modest effect. 
At the same magnitude limits but near the CCD center, the galaxy count correction is only 
^10%. Fortunately, our final galaxy cluster catalogs are insensitive to the application of these 
corrections (i.e., the same clusters are detected whether or not we apply these corrections). 



3.4. Astrometric Calibration 



A precise astrometric calibration is required to obtain accurate celestial coordinates for 
subsequent spectroscopic programs and to perform accurate catalog matching in order to 
generate V4—I4 colors. To perform this calibration, we first construct a global cartesian 
coordinate system for each color in each survey field by determining the relative offsets 
between the CCDs in the 4-Shooter as well as the offsets between adjacent scans. The global 
coordinate system provides the x-y position relative to the origin of the first scan in a survey 
field. The x-y offsets between the V4 band CCDs (CCDs f and 4) and between the I4 band 
CCDs (CCDs 2 and 3) are determined by using the photometric calibrations, or "fast" scans 
(see §3.2), which run perpendicular to the survey scans. The relative positions of objects 
on the fast scans give their true separations which, when combined with their apparent 
separations on CCDs f and 4 or CCDs 2 and 3, allow the inter-CCD offsets to be computed 
quite accurately. The inter-CCD offsets remain constant for any given scan but can vary 
from scan to scan as the 4-Shooter CCDs are periodically adjusted. Offsets between adjacent 
survey scans are determined by using the objects in the scan overlap regions. The precise x-y 
offsets between overlapping frames in adjacent scans also depends on the declination (since 
the scans converge as one moves towards the f 950.0 pole). To account for this, the observed 
offsets at the start, middle, and end of the adjacent scans are measured; the offset for each 
overlapping pair of frames is then calculated by a linear interpolation of the measured offsets. 
When the offsets are determined, the global x-y position for each object is computed and 
stored in the catalogs. Because the offsets are all derived using 4-Shooter data directly, 
the angular scale per unit position in the global system is identical to the CCD image 
scale. Indeed, we recover the published 4-Shooter CCD scale with great accuracy when our 
astrometric solution is computed. Our mean best fit 4-Shooter image scale (in arcseconds 
per pixel) is 0.33501 ± O.OOOII (^4) and 0.33503 ± 0.00021 (74). 

The astrometric reference objects are obtained by selecting 50-60 unsaturated [V4 ^ 
17) stars per survey field brighter than V4 = 19. The J2000 celestial positions of these stars 
are then measured from the digitized "Quick V" (epoch 1983) Palomar Observatory Schmidt 
plates (Lasker et al. 1990). The stellar centroid is measured from the digitized scan using 
the density-weighted position and the celestial coordinates are determined using the plate 
solution generated during construction of the HST Guide Star Catalog VI. (see Lasker et 
al. 1990 and Russell et al. 1990 for details). An astrometric solution for each 4-Shooter 
survey field is then computed using a 6-term per coordinate polynomial fitting procedure, 
with one iteration to exclude reference stars with large (> I''5) residuals. The typical rms 
deviations in the final solution are 0''48 in right ascension and 0''36 in declination. The 
scatter in the astrometric solution is due to two roughly equal components — uncertainties 
in measuring the relative scan and CCD offsets and uncertainties in the measuring the 
centroids of reference stars from the digitized Schmidt plates. 



4. CLUSTER CATALOG CONSTRUCTION 
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4.1. The Cluster Finding Algorithm 



The detection of a rich cluster at low redshift is successfully accomplished by many 
different methods {e.g., visual, percolation, wavelets, matched filters). For statistical studies 
of distant clusters, though, the important properties of the algorithm are the detection rate 
for small amplitude fiuctuations (e.^., poor or distant cluster completeness), the false positive 
detection rate (spurious detections due to chance projections on the sky), and the accuracy 
with which the selection function can be determined. The latter is particularly important 
because the stability and accuracy of the cluster selection function are direct measures of the 
objectivity of a given cluster detection procedure. Our approach to cluster detection is to 
make use of both positional and photometric information simultaneously by using a matched 
filter algorithm which optimally enhances the contrast of a cluster galaxy distribution with 
respect to the foreground and background galaxy distributions. Since the shape of such a 
filter will naturally depend on the redshift of the clusters being searched for, the amplitude 
of the filtered signal provides a useful redshift estimate. 

The matched filter method has many attractive features — 1) it is optimal for identifying 
weak signals in a noise dominated background (exactly the case for distant clusters), 2) 
photometric information is incorporated in an optimal way, 3) the contrast of overdensities 
which approximate the filter shape is greatly enhanced, 4) improved suppression of false 
detections, 5) redshift estimates are produced as a byproduct, and 6) good performance is 
seen over a large dynamic range in redshift. There are some disadvantages, in particular, 
one must assume a form for the cluster luminosity function and radial profile. We find that 
with sensible choices, however, the filters are broad enough to allow detection of systems 
with a range of properties. 

The objective is to design a filter which preferentially suppresses galaxy fiuctuations 
which are not due to real clusters. We can derive such a filter directly from a maximum 
likelihood estimator by modeling the spatial and luminosity distribution of the galaxies 
within a cluster and then determining the likelihood C that our model fits the actual data. 
We begin with the distribution of the number of galaxies per unit area per magnitude in any 
given patch of sky. This distribution is represented as 



counts (num/mag/arcsec^), P[r/rc) is the projected radial profile of the cluster galaxies 
(num/arcsec^), (f>[m — m*) is the differential cluster luminosity function (num/mag), and 
Ac/ is proportional to the total number of cluster galaxies and is a measure of the cluster's 
richness. The parameter m* is the apparent magnitude corresponding to the characteristic 
luminosity of the cluster galaxies and Tc is a characteristic cluster scale length (as seen in 
projection). 



_D(r, m) 



background + cluster 
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The likelihood C of the data given the model or, equivalently, the likelihood of having 
a cluster at a given position is given by 



In £ oc / In cr + 

JArea,m J Area 



[h{m) + AciP{r /rc)(f>{m — m*) — D{r, m)]" 
^2 



12) 



The above equation is valid when 6(m) is large enough that the Gaussian approximation 
holds; therefore, the statistics are dominated by the background galaxy counts and attribute 
all uncertainty to Poissonian noise in the background. This implies = b[m). We explicitly 
ignore the two point galaxy-galaxy correlation function. We do not solve for b[m) but assume 
it a priori. 



The likelihood is, thus, a function of 



m , r^, au' 



d A, 



/, I.e. 



jC 



jC(m*,rc,Aci) 



The 



dependence on m* and Vc allows one to obtain a redshift estimate for each cluster. We 
choose to accomplish this by maximizing the likelihood with respect to the characteristic 
m*. The characteristic magnitude, m*, is much more sensitive to the redshift than the 
cluster scale, Vc- 

Expanding the right hand side of the likelihood equation (Eq. 12), we get 
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The expression for Ac/, determined by setting ^ 
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/ P'^[r/rc) ^ d'^r dm 

Inserting A^ (Eq. 14) into the expression for the likelihood (Eq. 13), we find 

|/ -^(''/'''')'?^(™~™ )D{r,m) _ J P (^r / v c) (f>{m — m*) d^r dm 
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The expression in the numerator is simply the square of the excess above the background, 
as the background contribution is removed by the second term in Eq. 15. Since this term is 
a constant (if the integrals converge; see below), we could simply maximize the first term if 
we needed only to maximize the numerator. The denominator is not a constant; however, it 
is monotonic function which is proportional to some multiple of l/6(m*), which to a good 
approximation is just a power-law in the fiux. Thus, maximizing the expression on the RHS 
of Eq. 15 is roughly equivalent to maximizing 



. (p(m — ) 2 
In C oc I P[r/rc) — ^ D[r^m) d r dm 
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except for a shift Am* in the location of the maximum which is constant if b[m) can be 
accurately modeled as a power-law in m. We choose to use this expression not primarily 
because of its simpler form but because it depends only on the shape of the background, 
while Eq. 15 is very sensitive to the correct normalization of the background. In practice, 
there are no quantitative differences in the cluster detections (and their estimated properties) 
when maximizing the formal likelihood equation (Eq. 15), compared to maximizing Eq. 16. 

The apphcation of the matched filter to an input galaxy catalog is accomplished by 
evaluating the integral in Eq. 16 at every point in the survey and over a range of m* and Vc 
values (or equivalently, over a range of redshifts). The characteristic luminosity, Z*, and the 
intrinsic cluster scale length are assumed to remain fixed in physical units. The observables, 
m* and Tc, are assumed to vary with redshift as prescribed by a Eriedman- Robertson- Walker 
cosmology with Ho = 75 km s~^ Mpc~^ and qo = 0.5. The redshift dependence of m* also 
includes a k-correction discussed in §4.2. Because the galaxy distribution is a set of discrete 
coordinates and fiuxes, one can represent -D(r, m) as a series of (5-functions at the observed 
positions and magnitudes. The integral in Eq. 16 is, therefore, equivalent to the sum over 
all galaxies in a field 



The "filtered" galaxy catalog is created by evaluating this sum for each element in a two- 
dimensional array, S{i^j). The elements correspond to a series of grid points covering the 
entire survey area. The array, S'(z,j), can be treated as an "image" and we shall hereafter 
refer to its elements as pixels. In Eq. (17), Ng is the total number of galaxies in the catalog, 
P{rk{i, ])) is the angular weighting function (hereafter the "radial" filter), rk{i,i) is the 
separation between the k^^ galaxy and the (z,j) pixel center, nik is the magnitude of the 
ki\i galaxy, and L{mk) is the luminosity weighting function (hereafter the fiux filter). The 
dimensions of the array S{i^j) are dependent on both the angular size of the survey and the 
redshift to which the filters are tuned. The best choice of pixel size is to maintain a constant 
metric width that is approximately half a cluster core radius (tc), although the results are 
quite similar with pixels ranging from 0.25rc to 2rc. While the above summation runs over 
all galaxies, in practice the radial filter will have a cutoff radius and, thus, the execution 
time is greatly reduced for a large catalog. 

The optimal fiux filter is '^^7(m7 ^ ' however, the integral of this function diverges for 
Schechter (1978) luminosity functions with a < —1 as one sums to magnitudes fainter than 
m*. We can correct for this and, at the same time, place most of the weight of the filter 
near m* (where the contrast of the cluster against the background is the largest) by using 
the modified fiux filter 



cf){m — m*) 



D{r, m) cPr dm 
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Here we have introduced a power-law cutoff of the dimensionless form 10 ^^'"^ ™ ^ to keep 
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the function integrable. We choose /3 = 0.4 as this would represent weighing by the flux of 
the galaxy. 

We choose as the radial filter 

Pir Itc) = I = , if r < Tco 

otherwise (19) 

where is the cluster core radius and r^o is an arbitrary cutoff radius. The function P{r jvc) 
is circularly symmetric. However, since the diameter of the filter extends well beyond a 
core radius, systems which are not circularly symmetric are still easily detected, albeit at a 
different significance level. 

We choose the following normalizations for the radial and fiux filters 

P(r/rc) 27rr dr = 1 (20) 



<5{m-m*) dm= (^(m - m*) 10-°-^("-"') = 1 (21) 

Jo 

The limits of integration for P[r/rc) are effectively [O,rco]. The limits for L[m) are the 
magnitude range of the sample [O^niUm]- Figures 9 and 10 show the final radial and fiux 
filters function of filter redshift. 

With the above normalizations, the "background level" of a filtered image is S{i^j) = 
Sbg = 1. We expect the distribution of S{i^j) values for a given field (and a given filter 
redshift) should be a roughly Gaussian distribution with a standard deviation characterized 
by the Poissonian error in the background distribution, a mode of 1, and a tail toward high 
values of S{i^j) representing significant cluster detections. In practice, the mode (Sbg) is not 
strictly 1 but varies slightly with the filter redshift because we have estimated the background 
counts, 6(m), as the total galaxy counts in the field (background + clusters); we have no way 
of separating the two components a priori. Figure 11 shows the distribution of S{i^j) values 
for the 0^ field with the matched filter "tuned" to a redshift of z = 0.6. 



4.1.1. Effect of the Normalization on the Redshift Estimate 



The normalization of the fiux filter attempts to compensate for the portion of the fiux 
filter that extends beyond the magnitude limit of the sample. This correction is strictly 
appropriate for a pure background distribution h{m), as it assumes that the cluster galaxy 
counts are small compared to the field galaxy counts. However, when a significant section 
of the fiux filter is truncated by the survey magnitude limit, a small but correctable bias is 
introduced in the estimated redshifts derived from the matched filter for the most distant 
clusters. The redshift estimate for a cluster candidate is determined by finding the filter 
redshift at which the signal, S'(z,j), is maximized. The bias is essentially insensitive to the 
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radial filter normalization and we will ignore the radial filter P[r/rc) in the computations 
here. The fiux filter normalization then implies the following for the data D = b -\- (pel 



jm.^ dm _ f^'- ^^^^ dm 



+ 1 22 



This is explicitly correct for background only [(pel = 0), as S{i^j) = Sbg = 1. For real data 
{(pel > 0), the above normalization results in an overcorrection of the cluster signal at high 
filter redshifts; that is, 

JO b ^ -JO b /90\ 

jmn™ ^ joo ^ { ) 

This overcorrection becomes important only when a significant portion of the fiux filter 
extends beyond the magnitude limit of the sample (m* — )• mum)- 

We can quantify the effect by analytically calculating the cluster signal, S'c/, as a function 
of both the filter m* (or equivalently filter redshift Zest) and the cluster redshift. We calculate 
Sci from the expression 

^(m - m*;)$(m - m*) , /""^ n / a n/ / m2 

Sol = / — rr^ -dm / P4r)P{r/rc)d^r (24) 

Jo b[m) Jo 

where (p[m — m*i) and Pci{r) are the characteristic luminosity function and radial profile of 
the input cluster. We have assumed a cluster luminosity function which has the same form 
and k-correction as the fiux filter and a cluster radial profile which has the same form as 
the radial filter (see §4.2, §4.3, and Table 3 for explicit information on the filter parameters 
used). P[r/rc) and $(m — m*) are normalized according to Eq. (20) and (21). 

The dependence of Sd on filter redshift is affected by the assumed functional form of 
the background galaxy distribution b[m) and the magnitude limit of the sample. Figure 12 
shows Sci as a function m* (the equivalent filter redshift, Zest^ is listed along the top axis) for 
the V4 and I4 bands. Fach curve represents a cluster at a redshift from 0.2 to 1.2 in intervals 
of 0.1 (from top to bottom). The filter is run at each Zest (or m*), and the resulting values 
of Sci are plotted. The survey magnitude limits are indicated by dashed lines. 



For clusters at lower redshifts, Sd reaches a maximum at or close to the actual redshift 
and, hence, a reliable redshift estimate is obtained. However, for a cluster at 2; > 0.6, i.e., 
as m* — )• miira, the maximum of Sd cannot be used to gauge the redshift of the cluster 
since Sd is always maximized at the highest filter redshift. This behavior is the direct result 
of the overcorrection of the cluster signal by the adopted fiux filter normalization. When 
m* « miira, the problem is negligible. 

We can correct for this bias by applying an appropriate, redshift-dependent cluster 
signal correction [CSC) to the detected cluster signal, S{i^ j) — Sbg- We derive this correction 
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by making Eq. (23) into an equality, which is the case when the flux fllter is not truncated 
by the magnitude limit, i.e., 

CQ,C = ^ am / Jo j, am 

^rnum ^ dm / dm ^ ' 

The corrected signal is then 

' ~ CSC 

The correction factor is dependent on the actual cluster redshift and the assumed luminosity 
function of the cluster (pel (here we have assumed that the clusters have the same luminosity 
function and k-correction as the flux fllter), but not on the cluster richness. In deriving the 
CSC, we have assumed that the clusters have the same luminosity function, k-correction, 
and radial proflle as our matched fllter. The flducial redshift we have chosen for computing 
the CSC is 2; = 0.7. The CSC is well approximated by the expression (1 + z'^gfY ■ The 
advantage of this approximation is strictly for computational convenience. The cluster 
redshift estimates do not change signiflcantly if we use the full expression given in Eq. 
25. The best-flt functions are 

CSC{V,) = (1 + z^^f^ (27) 
CSCih) = (1 + 4?)^-° (28) 

Note that these corrections are only estimates to the cluster signal as they explicitly depend 
on the actual redshift of the cluster, its luminosity function, and its k-correction. Eigure 13 
shows Sci as a function m* (and Zest] top axis) after application of the CSC . The signal now 
peaks at or near the actual redshift (up io z = 1.2) and, thus, is a better redshift estimator 
(see §4.2.1 for a discussion of the redshift accuracy). Nonetheless, estimated redshifts of 
clusters at 2; ^ 0.7 in the band and z ^ 1.0 in the /4 band are still quite uncertain as the 
matched fllter is fltting only the bright end of the luminosity function. 



4.2. Cluster Detection Criteria 



The algorithm is run on each of the ten galaxy catalogs (flve survey flelds in two 
passbands), tuned to redshifts in the range 0.2 < z < 1.2 in 0.1 intervals. The required 
inputs are lists of galaxy positions and magnitudes. At each redshift a flltered image of 
the galaxy catalog is produced according to the prescription in the previous section. Use 
of a matched fllter does not remove the need to select a detection threshold. Detection 
is done on the flltered galaxy maps where the cluster contrast is improved. Clusters are 
identifled by searching for the local maxima within a moving box which is 1.25h~^ Mpc 
across. We center the box on each "pixel" of the flltered galaxy map array, S{i,j). A 
candidate cluster is registered when the central pixel in the box is the local maximum and it 
lies above a prescribed threshold. The detection threshold is taken to be the 95*^ percentile 
(approximately the 2cr level). 
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The resulting V4 and I4 cluster candidates are matched with each other to locate those 
systems which are detected in both bands. We then impose the additional criteria that the 
peak signal (highest valued pixel in a given detection) must exceed the 3(7 threshold in at 
least one passband, if the cluster candidate is detected in both passbands, or must exceed 
the 4(7 threshold if it is only detected in a single passband [note, we use the S'^f'' (Eq. 26) 
distribution to determine the peak signal and significance of all detections]. Detection in a 
single bandpass is often explained by the fact that the V4 scans cover a bit more area than 
the I4 scans or by the presence of a bright red or blue star which reduces the effective I4 
or V4 survey areas differently. There are, none the less, some interesting cluster candidates 
detected in only one of the two bands. These are listed or described in Tables 4 and 5. A 
cluster candidate must satisfy the above criteria in at least two of the matched filter redshift 
intervals to be included in our main catalog. This last criterion provides additional rejection 
of spurious detections. 

To compute the 95*^ percentile detection threshold and the significance of the peak signal 
in a cluster candidate, we examine the binned distribution of the values in the filtered galaxy 
map array, S{i^j). The distribution is typically Gaussian with a tail toward higher values, 
representing the cluster detections (see Figure ff). We characterize this distribution by 
calculating the background level, Sbg^ as the mode and the dispersion around the background, 
(76p, as 0.74f X [Q3 — Qi) where Qi and Q3 are the first and third quartiles, respectively. The 
coefficient is correct for computing the standard deviation from the interquartile range of a 
normal distribution. The significance of a given peak signal, Sp^ is then simply [Sp — Sbg) / (Jbg. 

The redshift estimate for a given cluster candidate is determined by finding the filter 
redshift at which the peak signal is maximized (e.^., see Figure f5). The choice of the k- 
correction and of the Schechter (f976) luminosity function parameters used to define L[m) 
strongly affect the redshift estimate but have little effect on cluster detection reliability 
providing extreme values are not used and that the luminosity functions of distant clusters 
are closely approximated by the Schechter form. We quantify this last statement in §4.3. 

4-2.1. Redshift Estimation Accuracy 

Figure f 4 shows the offset between the estimated and true redshift as a function of true 
redshift for ~ 400 simulated clusters. The luminosity functions of the simulated clusters 
exactly match the luminosity function used to compute the fiux filter function. The offset 
is small and is independent of the cluster redshift as long as the assumed k-correction and 
form of the cluster luminosity function are accurate and the C SC has been applied. 

The real test of the accuracy of our redshift estimates, of course, is to compare them 
with the observed redshifts of real clusters in our survey. Unfortunately, there are only 
six clusters in the survey which have reliably determined redshifts, although they span 
a reasonable range in redshift, 0.3 < Zohs < 0.8. Figures f5a and f5b show the offset, 
Zest — Zohs^ as functions of the matched filter signal {Sli '') and passband for these six clusters. 
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In all cases, the matched filter signal reaches its peak within Zobs i 0.2. A subregion of the 
filtered 13^ galaxy catalog (with the CSC applied) centered on the GHO cluster 1322+3028 
{zobs = 0.70) is shown in Figure 16 at the 11 filter redshift intervals (0.2 — )• 1.2) used. 
The GHO cluster 1322+3027 {zobs = 0.75) also lies in this subregion. The most distant 
clusters appear to have estimated redshifts which are systematically lower (by 0.1 — 0.2) 
than their observed redshifts. Such a trend is expected if we have underestimated the amount 
of spectral evolution (i.e., overestimated the k-correction). Since our nominal fiux filter is 
constructed assuming no evolution (see §4.3), this trend is not surprising. However, given 
that precise evolutionary corrections are not well known, we simply publish our redshift 
estimates uncorrected for such systematic trends and present this result as a caveat to users 
of this catalog. 

The V4 and I4 redshift estimates are quite consistent. Of the 74 matched clusters in 
our main catalog, 58 clusters (78%) have \zv — zi\ < 0.2. Of the 16 matched clusters 
with \zv — zi\ > 0.2, half are cases of distant clusters located near clusters at moderate 
redshifts where the detection algorithm locked onto the distant cluster in one passband and 
the foreground cluster in the other passband. The remaining differences are consistent with 
the inherent errors in the redshift estimates. 



4-2.2. Richness Measure 



Once we have determined the estimated redshift of the cluster, the cluster richness is 
given by calculating at Zest- From Fqs. 11 and 14 and the adopted filter normalizations 
(Fqs. 20 and 21), A^ is given by 



S{^.J)-l 

foo P2(r/r \ ^Trr rlr ["^l^m 4>^ (m-m*) f (m) j 
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where /(m) = 10"°-'*™. We use S{i^j) prior to application of the CSC (the CSC is only 
needed for redshift estimation). The normalization of the differential cluster luminosity 
function (f>[m — m*) is / cf){m — m*) f{m) dm = L* (which is exactly equivalent to Fq. 21); 
this normalization implies that A^ is a measure of the effective number of L* galaxies in a 
cluster. We can easily see this by determining the total fiux of the cluster (Ltot)- Integrating 
Fq. (11) over luminosity [L = 10"°''*™), we have 



Ltot = J 10~'^-^"' D{r,m) d\ dm 

= I 10-°-4" b{m) dm I d'^r + A^ j (p{m - m*) IQ-O-^" dm j P{r/r,) d'^r 

= constant + AdlO-°-^"' (30) 
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Therefore, the cluster luminosity is L^i = Ac/10 = AdL*, implying that Ad is the 

equivalent number of L* galaxies in the cluster. 

4-2.3. Ahell Richness Measure 

Since estimates of cluster richness based on parameters from the matched filter cluster 
algorithm, such as the peak signal or the detection area, depend on the cluster profile and 
estimated redshift, we make an independent estimate of richness from the cluster galaxy 
counts in each of the two passbands. We follow the specification of Abell (1958) and calculate 
richness (Nr) as the number of member galaxies (above background) that are brighter than 
ms + 2™ (where ms is the magnitude of the third brightest galaxy) and located within a 
projected radius of 1.0h~^ Mpc from the cluster center. We use a radius of 1.0h~^ Mpc, 
rather than Abell's historical 1.5h~^ Mpc, because of the growing uncertainty in background 
subtraction at larger radii. Even at 1.0h~^ Mpc, it is still difficult to calculate ms accurately 
from the cluster number counts; therefore, we first calculate the central density {Nq; Bahcall 
1981). That is, we determine ms from the background subtracted galaxy counts within the 
smaller radius of 0.25h~^ Mpc. 

For each field, we calculate the differential background galaxy surface density 
(galaxies/deg^/0.20 mag) by using regions which contain no detected clusters. For each 
cluster in the field, we bin all galaxies within the central 0.25h~^ Mpc (at the estimated 
cluster redshift) in the same 0.20 magnitude bins as the background galaxies. The expected 
number of background galaxies in each magnitude bin is subtracted. From the final cluster 
number counts, we determine the magnitude bin which contains the third brightest galaxy. 
The sum of all the galaxies in bins with m < ms + 2™ is defined as the central density of the 
cluster. We then perform the same procedure at a radius of 1.0h~^ Mpc, except that we fix 
ms at the value determined in the central density calculation. 

With the above method for estimating ms, we find that there is a good correlation 
between Nr computed within a radius of 1.0h~^ Mpc and Nr computed within a radius 
of l.bh-^ Mpc. Specifically, Nr{1.0)/Nr{1.5) = 0.72 ± 0.05. We use this relationship to 
estimate the conventional Abell richness of a distant cluster candidate in order to make 
some comparisons with the space density of low redshift clusters in §5. 

Figure 17 shows the relation between Nr[1.0) and for both the real cluster candidates 
and for cluster detected using the algorithm in simulated galaxy catalogs. The correlation 
between the two richness parameters is significant at a greater than 99.2% level (using the 
Spearman rank correlation function), indicating that these richness measures are meaningful. 
The rms scatter, however, is large with sigmaA^i ~ 23 for a given Nr. The scatter is, in part, 
due to the dependence of A^ on the shape of both the matched filter and the actual cluster, 
as well as background subtraction and redshift errors in the calculation of Nr. When ms + 2 
is fainter than the survey magnitude limit, the estimated cluster richness count is highly 
uncertain and is often too small. Consequently, for the highest redshift clusters (indicated 
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by stars in Figure 17), the richness estimate Nr systematically underestimated the actual 
cluster richness. 

4.3. The Cluster Catalog 

The matched filter and cluster detection parameters used to generate the cluster catalog 
are listed in Table 3. These parameters were chosen because they provide a reasonable 
compromise between minimizing spurious detections and maximizing completeness at high 
redshift. We choose 100h~^ kpc and lh~^ Mpc as the nominal values for the core radius and 
the cutofi' radius, respectively, in the definition of the radial filter (Eq. 19). This choice was 
motivated by several studies on the profiles of nearby [z < 0.2) Abell clusters (e.^.. Beers & 
Tonry 1986; Oegerle et al. 1987; Postman, Geller & Huchra 1988; Broadhurst 1995; Smail et 
al. 1995; Squires et al. 1995; Tyson & Fischer 1995). Quantitative descriptions of how the 
selection of radial filter parameters affects catalog membership are given in §4.4.2 and §4.4.3. 
A detailed analysis of the radial profiles of the distant clusters found in the present survey 
will be presented in paper II (Lubin & Postman 1995). Briefiy, we find that the typical 
surface density profile of our cluster candidates has Tc ^ 140/i~^ kpc and a power law profile 
on scales out to at least lh~^ Mpc with slope —1.4 it 0.4. 

The luminosity function parameters used are based on observations of low redshift 
clusters and field galaxies (Kirshner et al. 1979; Lugger 1986; Postman et al. 1988; CoUess 
1989; Marzke et al. 1994). We assume that the shape of the luminosity function is 
independent of redshift and adopt a k-correction appropriate for a non-evolving elliptical 
galaxy (see Figures 3a and 3b). While there is evidence to suggest that the galaxy luminosity 
function does depend on redshift (Fales 1993; Lilly 1993; Lonsdale & Chokshi 1993), our 
choice to ignore this primarily affects the accuracy of the redshift estimates (see §4.2.1) and 
does not significantly affect the catalog membership. This is because catalog membership 
is more sensitive to the shape of the fiux filter than the location of its peak. The redshift 
estimate, on the other hand, is more sensitive to the location of the peak of the filter. 
Plausible changes in the shape of the cluster luminosity function with redshift result in 
changes in the fiux filter shape which are significantly less than its width. For example, if we 
vary the slope of faint end of a Schechter luminosity function from a = —0.6 to a = —1.6 the 
FWHM of the fiux filter varies by ^ 15% relative to the FWHM of our nominal (a = —1.1) 
filter (the fiux filter becomes broader as we steepen the faint end slope). The peak of the fiux 
filter shifts from 0.2 mag brighter than the peak of our nominal filter to 0.3 mag fainter than 
the peak of the nominal filter as we vary a over the above range (the filter peaks at fainter 
magnitudes as we steepen the faint end slope). For comparison, the FWHM of our nominal 
fiux filter is 1.87 mag, which is significantly larger than the above shift. To be sure, we ran 
the algorithm on the 13^ field and varied the slope of the luminosity from —0.6 to —1.6 in 
0.2 intervals. The resulting catalogs did not vary by more than 10% in size (relative to the 
13^ catalog published here) and the detailed agreement is excellent (17 of the 19 13^ cluster 
candidates were found at all the a values explored). The width of the fiux filter is insensitive 
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to the characteristic luminosity, Z*, or the k-correction. The location of the peak of the flux 
filter, however, is quite sensitive to these parameters. If the amount of spectral evolution 
is underestimated (i.e., the k-correction is overestimated) then the redshift estimate will be 
underestimated and visa versa. For example, if z ^ 0.5 cluster galaxies have spectral energy 
distributions similar to non-evolving Sbc galaxies then using a non-evolving elliptical galaxy 
k-correction in the construction of the fiux filter will result in Zest — Zohs ~ —O.I at Zohs = 0.5 
and Zest - Zobs ~ -0.2 at Zobs = 0.8. 

Table 4 contains the 79 cluster candidates which satisfy the detection criteria with the 
above matched filter parameters. For each cluster, the first column contains a numeric ID, 
columns 2 and 3 contain the J2000 equatorial coordinates of the cluster (derived from the 
weighted average of the locations of the peak signal in the V4 and I4 bands for systems 
detected in both bands), columns 4 through 10 contain, respectively, an effective V4 band 
cluster radius in arcseconds, the V4 band redshift estimate, the peak V4 band cluster signal 
(S^i^)^ the significance of the peak signal, the richness estimate (Ac/) derived from the V4 
cluster signal, and two richness estimates derived from the unfiltered V4 band galaxy catalog 
[No is the background subtracted galaxy count within 0.25h~^ Mpc; Nr is the counts within 
1.0h~^ Mpc; see §4.2.3). Columns II through 17 contain the same information for the I4 
band. The effective radius is defined as 



where A^et is the area of the cluster detection that lies above the 95 percentile threshold. 
If a cluster is only detected in a single passband, we measure its properties in the undetected 
band using the position and estimated redshift derived from the actual detection. Such 
clusters are identified in Table 4 by a missing entry in the effective radius column for the 
band in which the cluster was not detected. There is one cluster, ^^028, for which I4 data 
does not exist. This cluster is a low z (zest = 0.2) system which would have been easily 
detected in the I4 band. 

Single band detections with peak signal < 4(7, Zest ^ 0.4, and which appear to be 
real clusters via visual inspection of the CCD images are listed in our supplemental catalog 
(Table 5). The supplemental catalog is not intended to be complete but contains 28 excellent 
distant cluster candidates, one of which is a known GHO cluster. The notes to Tables 4 and 5 
contain details about some of the cluster candidates in the main and supplemental catalogs, 
including alternate identifications when available. 

Pictures of all 107 cluster candidates are displayed in Plates 1-6. The images are 
of the central 134" X 134" regions; north is up and west is to the right. The cluster ID is 
given at the top of each picture (supplemental cluster IDs are followed by the letter S). The 
clusters are shown in the passband in which the significance of the peak signal is highest. 
Occasionally there are narrow gaps in the images - these are just the unusable areas between 
4-Shooter CCDs. 

Although the PDCS fields are not all completely interior to the original Gunn & Oke 
(1975) survey areas, there are 27 GHO clusters that were imaged by the present CCD survey. 




(31) 
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With the adopted matched filter and cluster detection parameters, 24 of the 27 GHO clusters 
are detected but only 16 of these satisfy all the criteria for inclusion into the main catalog. 
GHO cross identifications are given in Table 6. The 3 GHO clusters not detected at all are 
CL0030+0536, CL0030+0512, and CL1605+4106; they are all missed because their peak 
signals fall below the 95*^ percentile detection threshold. 

4.4. False Positive Rate 

The number of spurious clusters detected (i.e., the false positive rate) is a function of the 
detection threshold and the matched filter parameters. We are able, fortunately, to measure 
these dependencies with sufficient accuracy that quantitative results can be extracted from 
our cluster catalogs. In the following sections, we discuss how these factors affect cluster 
catalog membership. 

4.4.1. Peak Threshold Dependence 

To evaluate the dependence of the false positive rate on the detection threshold, we ran 
simulations on artificial fields using a simple Monte-Carlo technique to reproduce, as closely 
as possible, the projected galaxy and magnitude distribution of the actual PDCS fields. In 
order to reproduce our matching criteria in the cluster detection method, we created ten 
pairs of background galaxy fields by simulating each field as viewed through the V4 and I4 
passbands. We accomplished this by first generating a galaxy catalog in the V4 band by using 
all galaxies detected in the V4 band of one of the PDCS fields. The simulated field is divided 
into a grid of square cells of 6 arcmin sides. A background galaxy population is created by 
distributing these galaxies, with their corresponding magnitudes, at random within each cell 
with the observed variation in number from cell to cell and the average count per field. We 
then generate a matched simulated field in the I4 band by placing all the galaxies detected 
in the I4 band of the PDCS field at the positions of the galaxies in the simulated V4 field 
such that the galaxy at that position would have an appropriate color (see Figure 3b). 

The cluster detection algorithm is then run on each simulated background galaxy catalog 
using the same matched filter parameters as in the real data (see Table 3). We fix the 95*^ 
percentile detection threshold and standard deviation (cr) at values obtained from the actual 
galaxy fields. That is, we examine only those spurious detections in the uniform background 
which would have been considered cluster detections had they been present in the actual 
fields. Using the same detection criteria as discussed in §4.2, these spurious detections are 
first identified individually in V4 and I4 and then matched with each other to determine which 
are detected in both bands. The number of matched spurious detections (those detected in 
both passbands) as functions of peak signal threshold and radial extent (Eq. (31)) are shown 
in Figure 18. The histograms for the actual cluster detections are shown for comparison. 
There are 4.2 and 0.8 matched spurious detections deg~^ with peak signals of > 3(7 and 
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> 4(7, respectively. The observed density of matched cluster candidates is 13.4 deg~^ and 6.8 
deg~^ for peak signal thresholds of > 3(7 and > 4(7, respectively. Approximately 88% of the 
observed cluster candidates are matched detections. This would imply false positive rates of 
~ 31% and ~ 12% for the matched detections with peak signals > 3(7 and > 4(7, respectively. 
The observed densities of unmatched > 4(7 detections in the V4 and I4 passbands are 1.1 deg~^ 
and 0.6 deg~^, respectively. In the simulations, we find the densities of unmatched > 4(7 
detections in the V4 and I4 passbands are 0.7 deg~^ and 0.4 deg~^, respectively. This would 
imply a very high false positive rate for the unmatched detections (about 65%). Figure 18 
demonstrates, however, that the actual detections are typically larger than the spurious ones. 
If one restricts the radial extent to be at least 200 arcseconds, the false positive rates drop 
to ~ 15% and ~ 9% for matched detections with peak signals > 3(7 and > 4(7, respectively. 

We can use the Abell richness estimates (Nr) of our cluster candidates to check these 
estimated spurious detection rates. We do this by computing the percentage of clusters with 
richnesses which are within two standard deviations of zero (i.e. consistent with being simply 
background). If such clusters are all false positive detections then we find spurious detection 
rates of < 20% for > 3(7 matched detections, < 9% for > 4(7 matched detections, and < 22% 
for > 4(7 single band detections. These rates are in good agreement with our estimates of 
the false positive rates in the simulated field galaxy distributions discussed above. 

If the peak threshold is lowered below the 3(7 level, the number of detections increases 
significantly: the detection rate is about 50% higher if a 2.5(7 peak threshold is used. 
However, many of the additional detections are either superpositions of field galaxies or very 
poor systems. There are inevitably richer but quite distant clusters which fall below the 3(7 
threshold. Some are listed in the supplemental catalog (Table 5). Because the false positive 
rate at the 2.5(7 peak threshold is comparable with the probable number of real detections, 
we have adopted the more reasonable 3(7 peak threshold. In any case, the advantage of an 
objective detection algorithm is that a selection function can be accurately quantified and 
we do so in §4.7.1. 

4-4-2- Core Radius Dependence 

To quantify the dependence of the catalog size on the core radius used in the radial filter, 
we ran the cluster finding algorithm on the 0^ field using three different core radii: 100h~^ 
kpc (the nominal value), 33h~^ kpc (3 times smaller), and 300h~^ kpc (3 times larger). These 
tests were done with the cutoff radius held constant at lh~^ Mpc. All the candidates found 
when we use a radial filter with a 300h~^ kpc core radius are also found when the 100h~^ 
kpc core radius filter is used. Conversely, only 70% of the candidates found with the 100h~^ 
kpc core radius are found with the 300h~^ kpc core radius filter. Some of the candidates 
missed with the larger core radius filter are among our most distant cluster candidates. The 
number of cluster candidates found when we use a radial filter with a 33h~^ kpc core radius 
is twice the number found when a 100h~^ kpc core radius is used. As one lowers the core 
radius scale, one picks up smaller systems and consequently the spurious detection rate also 
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increases. Some of the smaller systems detected are poor groups of galaxies. The rest are 
almost certainly superpositions of field galaxies that are not physical associations. None of 
the systems detected with the 33h~^ kpc core radius filter but missed at larger radii are rich 
clusters. Simulations reveal as the adopted core radius is decreased the false positive rate 
increases, reaching nearly twice the rates quoted in the previous section when Vc = 33h~^ 
kpc. Conversely, the false positive rate declines to about 65% the rates quoted in the previous 
section when Vc = 300h~^ kpc but then so does the detection probability of compact systems 
and very distant systems. The selection of a 100h~^ kpc core radius thus appears to be 
appropriate for reliable identification of clusters out to 2; 1. 

4-4-3- Cutoff Radius Dependence 

The catalog size also depends on the cutoff radius of the radial filter. We have explored 
a range of cutoff radii from 700h~^ kpc to 3h~^ Mpc (with the core radius held constant at 
100h~^ kpc). The significance of a detection, on average, gradually decreases as one increases 
the cutoff radius from 0.7h~^ Mpc to 3h~^ Mpc, dropping by 20 — 40% over this range. The 
dependence of the detection significance on cutoff radius is less well defined, however, in the 
range 700h~^ kpc ^ ^ lh~^ Mpc. A cutoff radius of lh~^ Mpc appears to be about 
optimal and is indeed what we have chosen. The most robust catalog members are obviously 
the > 4(7 detections. Their inclusion in the final catalog is insensitive to the cutoff radius 
in the range explored here. The dependence of detection significance on the cutoff radius is 
understood. As one increases the cutoff radius, more weight is given to galaxies at larger radii 
and less to galaxies near the cluster center (a consequence of the radial filter normalization 
given in Eq. 20). If the cluster dimensions are significantly smaller than the cutoff radius 
then the matched filter signal will be reduced because of the increased contribution from 
interlopers at large radii. If the cluster dimensions are significantly larger than the cutoff 
radius, the signal may be truncated but this depends on the precise shape of the cluster 
profile (clusters with steeper profiles than the radial filter are not as significantly affected). 

4-4-4- Cluster Dehlending 

About 5% of the cluster candidates have a neighboring cluster candidate that lies within 
a projected separation of 2h~^ Mpc and that has an estimated redshift that differs by < 0.2. 
While it is possible (although unlikely) for clusters to be so close, an alternative explanation 
is that substructures within a cluster were detected as separate systems. For reference, fewer 
than 3% of the R > Abell clusters in the volume limited sample produced by Lauer & 
Postman (1994) have projected separations less than 2h~^ Mpc. In the notes to Tables 4 
and 5 we attempt to identify possible related systems. We have not explicitly removed such 
systems because it is possible that they really are separate clusters (visual inspection was 
inconclusive). Follow-up spectroscopic observations will eventually resolve such cases. In 
any case, the effect is significantly smaller than the estimated spurious detection rate. 
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4.5. False Negative Rate 

Real clusters that are not detected fall into two categories. First are those clusters or 
groups which fail to satisfy either the 95*^ percentile detection threshold or the peak signal 
threshold because they are either very poor systems or extremely distant systems. It is 
inevitable that we will miss such systems. Fortunately, we are able to precisely quantify this 
selection bias. Second are those clusters which fall below the detection thresholds because 
they exhibit very flat profiles. If such systems exist, they may be missed because their shape 
is not well matched to the shape of the radial filter. Poor clusters with very "cuspy" profiles 
(core radii ^ 30h~^ kpc) might also be missed as they will be significantly smoothed, although 
no presently observed clusters have such narrow profiles. An accurate determination of the 
cluster selection function is essential, therefore, not only for understanding the completeness 
of our survey but for interpreting the observed cluster abundances in the context of competing 
cluster formation models (e.^., Frenk et al. 1990). 

4.5.1. The Selection Function 

The probability that a cluster will be included in the sample depends on its redshift, 
richness, profile, the evolutionary history of its galaxies, and the observation passband. In 
order to derive the selection function which describes all these dependencies, we first need 
to establish the relationship between the matched filter signal and these quantities. To do 
this accurately, we simulated a total of 10560 clusters at different redshifts (0.2 < < 1.2), 
richnesses (Abell R=l— 4), projected radial profiles (from fiat to r~^), passbands [V4 and 74), 
and assumed two different k-corrections: one appropriate for a non-evolving elliptical galaxy 
and one appropriate for a non-evolving Sbc galaxy (see Figures 3a and 3b). The simulated 
clusters all have Schechter luminosity functions with faint end slopes that are identical to 
that used in the construction of the fiux filter (i.e., a = —1.1). Fach simulated cluster was 
then added into each real galaxy catalog and the cluster detection algorithm was run. This 
allows us to directly compute the matched filter signal and detection frequency as functions 
of redshift, richness, profile slope, k-correction, and passband. 

Figure 19 shows the matched filter richness, Ac/, as a function of redshift, richness, k- 
correction, and passband for simulated clusters with r"^'"* profiles (which is the average profile 
slope of the actual cluster candidates). The matched filter parameters used are identical to 
those listed in Table 3. This is an important calibration since the peak signal and are 
not related to the cluster galaxy overdensity in a simple way. The observed values of 
all the cluster candidates (including the supplemental clusters) are plotted for comparison, 
along with the typical detection threshold level. The richness of detected clusters tends to 
increase with redshift but as Figure 19 demonstrates the precise trend depends on the degree 
of spectral evolution of the cluster galaxies. 

We also use our simulations to quantify the effect of mismatch between the cluster profile 
and the radial filter. The range of cluster profiles used in the simulations (r° to r~^) bracket 
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the range of profiles typically seen in Abell clusters. It was found that, as expected, the 
clusters yield the strongest signal and the flat (r°) clusters yield the weakest signal. The 
relative matched filter signal ratios for r°, r~^, and clusters of a given richness class are 
0.20 : f : 2.85. These ratios are independent of richness class and redshift. The uncertainties 
in these ratios are about ±10%. In other words, a richness class I cluster with a fiat profile 
(if such a system exists) would yield a peak signal that is about 14 times smaller than that 
for a richness class I cluster with a profile at the same redshift. 

The cluster selection function is computed by simply counting the number of simulated 
clusters that are detected as a function of redshift, richness, profile, k-correction, and the 
observation passband. The selection functions for clusters with r"^'"* profiles are shown in 
Figure 20. The trends are well understood. Obviously, the richer clusters are detected more 
frequently than poorer ones at high redshift. Likewise, clusters with the Sbc k-correction 
are more easily detected at high redshift than those with the E-galaxy k-correction. Distant 
clusters are also more easily detected in the I4 band. The selection functions for clusters 
with fiat or profiles are qualitatively similar but, of course, differ quantitatively. 

4.5.2. Merged Clusters 

A standard problem with any two-dimensional or quasi two-dimensional cluster 
detection algorithm is separating 2 or more superposed clusters at different redshifts. Our 
simulations have shown that in only 10% of the cases where two clusters had a projected 
separation of 1.25h~^ Mpc or less, did the algorithm successfully detect both systems. This 
is largely a consequence of the 1.25h~^ Mpc detection box used to locate local maxima in 
the filtered galaxy maps. All clusters with projected separations greater than 1.25h~^ Mpc 
are resolved. By selecting a smaller detection box one will resolve finer scale structures 
but at the cost of introducing many more multiply detected clusters (see §4.4.4). Indeed, 
if we reduce the size of the detection box by 50% the number of multiply detected clusters 
increases almost tenfold. The frequency of occurance of superposed systems is a function 
of the spatial distribution of distant clusters. Since this is not yet known, it is difficult 
to estimate how many systems are missed or have richness overestimated because of this 
effect. Spectroscopic follow-up surveys of these systems are required to definitively answer 
this question. 

5. DISCUSSION 

One of the prime motivations for conducting this survey is to constrain the comoving 
density of clusters as a function of redshift. This will be accurately accomplished only when 
redshifts exist for a substantial number of these distant clusters. At the present time we 
have spectroscopic redshifts for only a handful and, thus, we present preliminary constraints 
based upon our estimated redshift data. 
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In Figure 21 we plot the cumulative number of clusters per square deg function 
of redshift and richness [Ad and NR[1.0h~^ Mpc)) for the PDCS. The A^(< z) for the 
GHO catalog, the Abell/ACO catalog, the APM cluster catalog (Dalton et al. 1992), and 
the Edinburgh-Durham-Milano (EDM) cluster redshift catalog (Nichol et al. 1992) are also 
shown. The Abell cluster curves are extrapolations normalized to the z < 0.07, R > Abell 
cluster space density and assume this density remains constant in comoving coordinates. For 
reference, a richness class 1 cluster should have 35 < Nr < 56. If the same cluster has z ^ 0.7 
and a r"^'"* profile, then it would have 25 ^ ^ 60 in the V4 band and 30 ^ A^ ^ 65 in 
the I4 band (see Figure 19). No corrections have been made for false positive detections in 
Figure 21. 

Figure 22 shows the observed and expected differential redshift histograms for the 
clusters in our main catalog with richness class > 1. Figure 23 shows the cumulative 
histograms. The observed histograms are based solely on our estimated redshifts and we have 
corrected the counts for false positive detections using the rates estimated from simulations 
(see §4.4.1). The expected distributions, for a given passband, are computed from the 
expression 

N,,p{z, k{z)) = V{z) J2{p,{z)P,{z, k{z))) (32) 

8 = 

where Nexp{z^ k[z)) is the expected number of clusters per square degree per unit redshift for 
an assumed k-correction model specified by k[z)^ ^{^) the volume of a comoving element 
from z to z -\- 6z subtending a 1 square degree solid angle, Pi{z) is the space density of 
clusters in this volume with richness class z, and Pi{z^ the probability that a cluster 

of richness class i at redshift z would be detected by our cluster selection procedure for a 
given passband. Since the selection probability is a function of projected cluster profile, we 
need to assume something about the profiles of clusters in order to make our prediction. 
The projected profile of a typical z < 0.2 cluster is at least as steep as but fiatter 
than over the range 100 ^ r ^ 1000/i~^ kpc (Beers & Tonry 1986; Oegerle et al. 1987; 
Postman, Geller & Huchra 1988; Broadhurst 1995; Lubin & Postman 1995; Smail et al. 1995; 
Squires et al. 1995; Tyson & Fischer 1995). For this computation, we adopt the selection 
function for a r"^'"* profile — the mean profile shape for the clusters in this catalog. Lastly, 
we assume a constant comoving density of clusters and a constant cluster mass function 
(i.e., the relative abundances of clusters of different richness classes remains independent 
of redshift). The actual density of the various richness classes is derived from the cluster 
luminosity function given by Bahcall (1979), normalized to match the observed space density 
of low redshift clusters. We compute the predicted cluster counts for two different galaxy 
population scenarios - one where clusters are made up of nonevolving elliptical galaxies and 
one where clusters are made up of nonevolving Sbc galaxies (see Figure 3). The different 
scenarios account for the two substantially different expectations as a function of redshift 
and passband seen the Figures 22 and 23, although in the I4 band cosmological effects are 
comparable with the effect of varying the k-correction. The estimated redshift distribution of 
the Palomar survey clusters is consistent with the hypothesis that the typical distant cluster 
galaxy is bluer than a nonevolving elliptical galaxy at the cluster redshift. 
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The predicted cluster counts depend critically on which data are used to determine the 
low redshift normalization. If we normalize to the observed space density of R > 1 Abell 
clusters, we find that the expected number of such clusters with z < 0.3 per square degree is 
about a factor of 5 (i2) smaller than the number in our Palomar catalog even after applying 
a statistical correction for spurious detections. This "discrepancy" is clearly demonstrated 
in Figure 21 — PDCS clusters with Nr > 35 or with ^ 40 should correspond to Abell 
richness class 1 systems. The PDCS certainly has better completeness for poor clusters 
than either the Abell or GHO catalogs although assessing the exact incompleteness of these 
catalog is difficult given the visual selection procedures used in their construction. The space 
density of clusters in the PDCS is consistent with that previously found by GHO; however, a 
direct comparison is not very meaningful since there is no richness information for the GHO 
clusters and the GHO cluster selection function cannot be easily quantified. Extrapolating 
the density of 2; ^ 0.1, R > 1 Abell clusters to 2; > 0.3 results in a substantial underestimate 
of the cluster space density inferred from the present survey. The uncertainty in the density 
of 2; < 0.3 clusters from the PDCS is large but improves significantly hy z = 0.6 where R 
> 1 PDCS clusters still outnumber the predicted R > 1 Abell/ACO cluster density by a 
factor of 5 (after false positive correction). Hence, the expected counts shown in Figures 22 
and 23 are based on a comoving density of R > 1 clusters that is 5 times larger than that 
derived from the Abell and ACO catalogs. Our assumed Abell cluster space densities are 
1.13 X IQ-^h^ Mpc-3, 4.04 x IQ-'^h^ Mpc'^, 1.28 x IQ-'^h^ Mpc'^, and 2.74 x lO'^h^ Mpc'^ 
for richness classes 0, 1, 2, and > 3, respectively. These densities are derived from Bahcall 
(1979) but are in good agreement with subsequent Abell cluster redshift surveys (Huchra 
et al. 1990; Postman et al. 1992; Lauer & Postman 1994) and southern hemisphere redshift 
surveys of systems comparable to R > Abell clusters (Dalton et al. 1992; Nichol et al. 
1992). 

At 2; < 0.6, the slope of the PDCS A^(< z) vs z relation is consistent with a constant 
comoving density, albeit a density 5 times higher than the R > 1 Abell cluster density. 
The fiattening of the slope of the A^(< z) vs z relation at 2; ^ 0.6 is primarily due to the 
PDCS selection function (as demonstrated by Figure 23). Figure 22 shows that a significant 
constraint on the spectral evolution of cluster galaxies is the modal redshift of the cluster 
sample — for clusters dominated by non-evolving elliptical galaxies, the predicted mode of 
the PDCS cluster redshift distribution is between 0.3 — 0.45 whereas for clusters dominated by 
evolving elliptical galaxies (or non-evolving Sbc galaxies), the predicted mode may be as high 
as 0.75. The uncertainties in the redshift estimates, however, prevent us from constraining 
the cluster space density at 2; ~ 1. Follow-up spectroscopic surveys should allow a reasonable 
determination of the behavior of the space density of clusters at 2; > 0.6. Postman et al. 
(1995) have completed a contiguous 13 square degree / band survey which is approximately 
1.5 mag deeper than the PDCS in an additional effort to address this issue. 

There are several possible explanations for the difference between the abundances of 
Abell and PDCS clusters. The Abell catalog could easily be missing half the clusters 
(presumably mostly R < 1). The objectively derived APM and FDM cluster catalogs yield 
a somewhat higher space density of clusters than that derived from the Abell catalog when 
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cut at comparable richnesses although it is not inconsistent with the Abell value. However, 
nearly 60% of the EDM clusters with mio < 18.5 are not detected in the AGO catalog. 
This is true at all richness cuts (Lumsden et al. 1992). Henry et al. (1992) produced an 
x-ray selected cluster sample to study cluster evolution. At 2; < 0.08, they detect 7 clusters 
but only 3 of these are in the Abell catalog. The 4 x-ray selected clusters not in the Abell 
catalog have x-ray luminosities that are comparable with the ones that are. X-ray luminosity 
is highly correlated with richness and central density (Edge & Stewart 1991). 

Previous observations have shown that clusters at 2; ~ 0.4 have galaxy populations 
which are significantly bluer than their low redshift counterparts (Butcher & Oemler 1984; 
Oemler 1992). At higher redshifts, Aragon-Salamanca et al. (1993) report that "...by z ~ 0.9 
there are no cluster galaxies as red as present day ellipticals." Such evolution would enhance 
the detection rate of high redshift clusters at a given threshold (particularly in the band). 
Recent HST image data (Dressier et al. 1994; Couch et al. 1994; Ellis 1994) reveal that 
many of the blue (</ — r < 1.2) galaxies in 2; ^ 0.5 clusters are either "normal" spirals or have 
peculiar morphologies, producing non-elliptical fractions which are 3 to 5 times higher than 
the average current epoch cluster. A number of physical mechanisms have clearly reprocessed 
the morphological mix between epochs. The dense regions of clusters do provide unique sites 
where processes like ram-pressure stripping and galaxy-galaxy collisions can, in principle, 
become efficient means of altering the properties of a significant fraction of the galaxy 
distribution. One explanation could be spiral destruction (Dressier 1980; Postman & Geller 
1984; Whitmore, Gilmore, & Jones 1993), although observations at redshifts approaching 
unity suggest that the fraction of excessively blue galaxies may jump to as high as 80% 
(Rakos & Schombert 1995). Eading through cessation of star formation may also play a role, 
but studies at low redshift indicate that this process alone cannot dominate the evolution of 
blue and/or spiral galaxies. Gunn & Dressier (^1988) find that the spectra of cluster galaxies 
with z ^ 0.6 show, on average, smaller 4000 A decrements and a higher frequency of post- 
starburst features (the "E-|-A" spectral class) than those at 2; < 0.6. It is likely that in 
the cluster environment a combination of starbursts, tidal destruction, fading, and stripping 
conspire to produce the observed time variation of photometric and spectroscopic properties. 

Some optically selected high redshift "clusters" may be larger scale structures 
{e.g., "Great Wall"-like features) viewed at orientations that make them mimic clusters. 
Contamination by large scale structures would increase both the number of detections and 
the blue galaxy component (sheets of galaxies at low z are composed largely of nonelliptical 
galaxies). Only an "edge-on" sheet would potentially be detected as cluster candidate since 
the galaxy surface density in a face-on sheet is too low (de Lapparent et al. 1986, 1989). 
If the fraction of sheets with near "edge-on" orientation is independent of direction or 
redshift then the probability of viewing a sheet of galaxies along its narrow dimension is 
just proportional to the volume of space surveyed. Hence, deep optical cluster surveys may 
be more contaminated by such structures than low z cluster surveys. Lefevre et al. (1994) 
spectroscopically verified the existence of a galaxy structure of intermediate scale (~ 2 — 3/i~^ 
Mpc) &i z = 0.985. 
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Our estimated spurious detection rate due to alignments of field galaxies is ^ 30% 
(§4.4.f), and we are confident that field galaxy superposition is not the primary source of 
the factor of 5 discrepancy in cluster counts. Couch et al. (f 99f ) claim that the higher cluster 
density in their J band survey could, in part, be explained by foreground contamination from 
spiral rich groups. However, 88% of our clusters are detected in both the and /4 passbands 
which argues against severe contamination from such groups. A related effect, which is 
probably important for our survey, is an enhanced detection rate due to the superposition 
of poor groups and/or clusters. Systems with redshifts that differ by A2; ^ 0.2 and with 
projected separations of ^ lh~^ Mpc would be detected as a single cluster with richness 
significantly higher than any of the individual superposed components. The Abell catalog 
contains a number of such cases {e.g., A4f9, Af3f8, Af63f). Spectroscopic follow-up is the 
only unambiguous way to measure the magnitude of this problem. We note, however, that 
the matched filter algorithm is able to resolve overlapping systems if the projected separation 
is > l.2bh-^ Mpc. 

Finally, we address the question of how our measuring errors might contribute to this 
discrepancy. The space density of clusters in the Abell catalog is roughly exponential with 
Abell's richness parameter. If / is the fraction of clusters at (or above) some richness 
i?, din f / dR — f.f; that is, the space density of clusters falls by a factor of ~ 3 for 
each successive richness class. The richness class is roughly logarithmic with total galaxy 
population (which is proportional to Ac/), i.e. din / dR 0.45. Therefore, the fraction of 
clusters above a given value of goes as a power law function of A^. 

f oc A;,'-^ (33) 

Thus, if we have overestimated the richnesses of all our clusters by a factor of ~ 2, 
the discrepancy can be explained. We have examined this possibility through extensive 
simulations of clusters with a broad range of profile shapes and color evolution (Lubin 
1995). Firstly, we examine the relation between the actual richness and Nr. For most 
of the very distant clusters, we severely underestimate their richness; however, for some of 
the nearby clusters {zest ^ 0.4), we can misclassify R = clusters as R = I clusters and 
vice versa due to uncertainties in background subtraction and redshift errors. In the worse 
case, these misclassifications could cause a factor of ~ 2 overestimate in the space density 
of clusters at 2; < 0.4. Secondly, we have examined the effect of the profile shape on the 
filter richness parameter A^. The magnitude of A^ is proportional to the steepness of the 
cluster profile. Complete resolution of the discrepancy would require Abell R > I clusters 
to have A^ > 85 and A^ > 70 in the V4 and I4 bands, respectively. These values would 
correspond to an average cluster surface density profile of or steeper. This profile is not 
consistent with the PDCS clusters. The average profile shape is r"^'"* (Lubin & Postman 
1995). Lastly, we have examined the effect of the shape of the cluster luminosity function on 
the relation between A^ and the Abell richness estimator. The relationship depends on the 
faint end slope, a, but not on the k-correction. For values in the range —0.6 > a > — f .3, the 
relationship varies by only ±30%. The trend is such that if the cluster LF slope is greater 
(less) than the nominal fiux filter value, a R = I cluster would yield a lower (higher) A^ 
value. 



A more subtle question is the effect of systematic or random redsliift errors. One miglit 
naively think that the error could be explained by a relatively small systematic underestimate 
of redshift as the volume, and hence the number, increases dramatically with distance. In 
the Euclidean case, dlnV/ din z = 3.0 (where V is the comoving volume), so a factor of 
five in space density would require a systematic underestimate of the distance by a factor of 
0.6. At a redshift of 0.5 in an = 1 Einstein-deSitter universe, the value of dlnV/ din z is 
smaller, about 2.2; to explain the discrepancy, we would have to underestimate the distances 
by slightly more than a factor of 2. This seems hardly possible; however, the situation may 
be more complex as our estimate of richness is also distance-dependent. If we underestimate 
the distance, we may underestimate the richness as we have assumed a brighter characteristic 
m*. The richness measured relative to this m* would be smaller than that measured relative 
to the actual, fainter m*. The result is that we would confuse very rich clusters at great 
distance, of which there are few, with poorer ones nearby. We can estimate the magnitude 
of the effect as follows. Erom the simulations, we find that in the vicinity of the correct 
redshift, the inferred richness changes with the assumed redshift like din Ad/ din Za ~ 2.1. 
Now let us focus on some assumed redshift Za and richness A"; and suppose that the clusters 
that we think are at that redshift are really at some other (say larger) redshift Zf. Let us 
then ask how the total measured number of clusters at z < Za and > A"; changes as the 
real redshift Zf changes. The number of clusters that we infer to be in our sample as we 
change Zf is the product of the total number of clusters, which we assume is proportional to 
the comoving volume, and the fraction of clusters at Zf which are sufficiently rich to have 
measured richness A"; at the inferred redshift Za- Thus 

dlnN dlnV dlnf 

CX "l" 

dlnz din z din z 

dlnV dlnf dlnA^i 

CX "l" 

dlnz J In Ac/ dlnz 
^ 2.2 - (2.4)(2.1) 
?^ -2.8 

This richness effect is opposite in sign to the volume effect and more than cancels the 
former; the numbers change with almost the inverse third power of the redshift. We need 
to have overestimated the redshift by a factor of about 1.8 to explain the discrepancy. This 
is extremely unlikely. The validity of this richness correction calculation depends on the 
origin of the redshift errors. As we have discussed in Sect. 4.3 and above, any systematic 
redshift error is most likely due to our assumption of no evolution, rather than the choice of 
an inappropriate luminosity function. Eirst, there is no evidence of gross variations in the 
luminosity function among clusters; second, our simulations, which are all based on a single 
Schechter luminosity function, display errors similar to the real sample and no substantial 
redshift biases. Therefore, we see no evidence to suggest that there are large luminosity 
function variations which would invalidate the notion of "richness" as measured by either 
the matched filter (Ac/) or the Abell technique (Nr) and the calculation for the richness 
correction (Eq. 34). 



(34) 
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Let us now consider the ("Malmquist") bias due to random errors. As long as the 
number-distance relation is concave up (which is the case for both the richness-uncorrected 
and the richness-corrected relation in Eq. 34), the effect of random errors in the distances 
is to cause one to overestimate the galaxy density in a cluster. If the distribution of errors 
in Inz is assumed to be Gaussian with standard deviation cr, the density overestimate is the 
factor 



Here 7 = din N / din z —2.8, so the factor is approximately F = exp(3.9(T^), which is 5 for 
cr 0.64 and 2 for a 0.42. At our mean cluster redshift of ~ 0.6, the former corresponds 
to a mean redshift error of ~ 0.4, the latter to an error of ~ 0.25. We believe the errors in 
our estimated redshifts to be ~ 0.2 (see section 4.2.1). Therefore, perhaps one-half of the 
discrepancy can be accounted for by the Malmquist effect, and it seems likely that evolution 
and incompleteness in the Abell sample can account for the rest. However, we must wait 
for a substantial number of real redshifts for this sample before it is possible to clarify 
this issue satisfactorily. We have already begun a number of related spectroscopic surveys. 
Such programs are time consuming and one goal of publishing this catalog is to provide the 
community with a starting point for further exploration. Over the next five years, results 
from these and other ground-based spectroscopic surveys plus imaging with HST and x-ray 
satellites should yield a substantial improvement in our understanding of the formation and 
evolution of galaxy clusters. 

These observations were made possible by the efforts of R. Lucinio, B. Zimmerman, 
M. Carr, E. Danielson, E. Lorenz, V. Nenow, J.D. Smith, and J. Westphal, who helped 
design, build, and test the 4-Shooter, and J. Carrasco, D. Tennant, M. Doyle, and J. 
Henning, members of the Palomar staff. We thank Jill Knapp for important discussions 
regarding star/galaxy classification, Robert Lupton for his excellent advice in statistical 
matters concerning the matched filter, and Tony Tyson for providing information about the 
response function of the Bj filter. We thank Bob Nichol, Michael Strauss, and the referee 
for their thoughtful and constructive comments on this manuscript. Robert Deverill and 
David Saxe created some of the software used for the data processing and analysis. This 



work was supported in part by NASA contracts NAGW-2I66 (MP), NAG5-I6I8 (DPS), 
NAS7-9I8 (JGH), NGT-5I295 (LML) and by NSE grant AST86-I8257A02 (JEG). 




d 



(35) 



- 36 - 



REFERENCES 

Abell, G. 0. 1958, ApJS, 3, 211. 

Abell, G. 0., Corwin, H. G., & Olowin, R. P. 1989, ApJS, 70, 1. 

Aragon-Salamanca, A., Ellis, R. S., Couch, W. J., & Carter, D. 1993, MNRAS, 262, 764. 
Bahcall, N. 1979, ApJ, 232, 689. 

. 1981, ApJ, 247, 787. 

Bahcall, N., & Soneira, R. M. 1983, ApJ, 270, 20. 
Beers, T. C, & Tonry, J. L. 1988, ApJ, 300, 557. 
Broadhurst, T. 1995, private communication. 
Burstein, D., & Heiles, C. 1982, AJ, 87, 1165. 
Butcher, H., & Oemler, A. 1984, ApJ, 285, 426. 

Coleman, G. D., Wu, C. C, & Weedman, D. W. 1980, ApJS, 43, 393. 
CoUess, M. 1989, MNRAS, 237, 799. 

Couch, W.J., EUis, R.S., Malin, D.F., & MacLaren, I. 1991, MNRAS, 249, 606. 
Couch, W.J., EUis, R.S., Sharpies, R.M., Smail, I. 1994, ApJ, 430, 121. 
Cousins, A. W. J. 1976, Mem.RAS, 81, 25. 

Dalton, G. B., Efstathiou, G., Maddox, S. J., & Sutherland, W. J. 1992, ApJ, 390, LI. 
. 1994, MNRAS, 269, 151. 

de Lapparent, V., Geller, M.J., & Huchra, J. P. 1986, ApJ, 302, LI. 
. 1989, ApJ, 343, 1. 

de Vaucouleurs, G. 1948, Ann.Astrophys., 11, 247. 
Dressier, A. 1980, ApJ, 236, 351. 

Dressier, A., Oemler, A., Butcher, H., & Gunn, J. E. 1994, ApJ, 430, 107. 
Eales, S. 1993, ApJ, 404, 51. 

Edge, A. C. & Stewart, G. C. 1991, MNRAS, 252, 428. 

Ellis, R.S. 1994, STScI Workshop on Quantifying Galaxy Morphology at High Redshifts. 

Erenk, C. S., White, S. D. M., Efstathiou, G., & Davis, M. 1990, ApJ, 351, 10. 

Griffiths, R. 1990, Hubble Space Telescope Wide Field & Planetary Camera Instrument 
Handbook Version 2.1 (Baltimore, Space Telescope Science Institute). 

GuUixson, C. A., Boeshaar, P. C, Tyson, J. A., & Seitzer, P. 1994, ApJS, in press. 

Gunn, J. E., Carr, M. L., Danielson, G. E., Lorenz, E. 0., Lucinio, R., Nenow, V. E., Smith, 
D. J., Westphal, J. A., Schneider, D. P., & Zimmerman, B. A. 1987, Opt. Eng., 26, 
779. 



- 37 - 



Gunn, J. E., & Dressier, A. 1988, in Towards Understanding Galaxies at Large Redshift^ R. 
Kron and A. Renzini, eds., Dordrecht, p. 29. 

Gunn, J. E., Hoessel, J. G., & Oke, J. B. 1986, ApJ, 306, 30. 

Gunn, J. E., & Oke, J. B. 1975, ApJ, 195, 255. 

Gunn, J. E., & Stryker, L. L. 1983, ApJS, 52, 121. 

Henry, J. P., Gioia, I. M., Maccacaro, T., Morris, S. L., Stocke, J. T., & Wolter, A. 1992, 
ApJ, 386, 408. 

Hoessel, J. G., Gunn, J. E., & Tliuan, T. X. 1980, ApJ, 241, 486. 
Jarvis, J. E., & Tyson, J. A. 1981, AJ, 86, 476. 
Kent, S. M. 1985, PASP, 97, 165. 

Kirshner, R., Oemler, B., & Scliecliter, P. 1979, AJ, 84, 951. 

Lasker, B. M., Sturcli, C. R., McLean, B. J., Russell, J. L., Jenkner, H., & Sliara, M. M. 
1990, AJ, 99, 2016. 

Lauer, T. R., & Postman, M. 1994, ApJ, 425, 418. 

Lefevre, 0., Crampton, D., Hammer, E., Lilly, S. J., & Tresse, L. 1994, ApJ, 423, L89. 

Lilly, S. J. 1993, ApJ, 411, 501. 

Lonsdale, C. J., & Cliokslii, A. 1993, AJ, 105, 1333. 

Lubin, L. M. 1995, Ph.D. Thesis, Princeton University. 

Lubin, L. M., & Postman, M. 1995, in preparation. 

Lugger, P. M. 1986, ApJ, 303, 535. 

Lumsden, S. L., Nichol, R. C., Collins, C. A., & Guzzo, L. 1992, MNRAS, 258, 1. 
Marzke, R., Huchra, J. P., & Geller, M. J. 1994, ApJ, 428, 43. 
Neushaefer, L. W., Windhorst, R. A., & Dressier, A. 1991, ApJ, 382, 32. 
Nichol, R. C., Collins, C. A., Guzzo, L., Lumsden, S. L. 1992, MNRAS, 255, 21P. 
Oegerle, W. R., Hoessel, J. G., & Jewison, M. S. 1987, AJ, 93, 519. 
Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713. 

Oemler, A. 1992, in Clusters and Superclusters of Galaxies^ A.C. Eabian, ed., Kluwer 
Academic Publishers, London. 

Postman, M., & Geller, M. J. 1984, ApJ, 281, 95. 

Postman, M., Geller, M. J., & Huchra, J. P. 1988, AJ, 95, 267. 

Postman, M., Huchra, J. P., & Geller, M. J. 1992, ApJ, 384, 404. 

Postman, M., Oegerle, W. R., Hoessel, J. G., & Lauer, T. R. 1995, in progress. 

Rakos, K. D., & Schombert, J. M. 1995, ApJ, 439, 47 



- 38 - 



Russell, J. L., Lasker, B. M., McLean, B. J., Sturcli, C. R., & Jenkner, H. 1990, AJ, 99, 
2059. 

Schechter, P. 1976, ApJ, 203, 297. 

Schneider, D. P., Gunn, J. E., & Hoessel, J. G. 1983, ApJ, 264, 334. 
Schneider, D. P., Schmidt, M., & Gunn, J. E. 1994, AJ, 107, 1245. 
Shectman, S. 1985, ApJS, 57, 77. 

Small, I., EUis, R., Eitchett, M. & Edge, A.C. 1995, MNRAS, 273, 277. 

Squires, G., Kaiser, N., Babul, A., Eahlman, G., Woods, D., Neumann, D.M. & Bohringer, 
H. 1995, ApJ, submitted. 

Thuan, T. X., & Gunn, J. E. 1976, PASP, 88, 543. 

Tyson, J. A. & Eischer, P. 1995, ApJ, 446, L55. 

Valdes, E. 1982, Proc. SPIE, 331, 465. 

Wade, R. A., Hoessel, J. G., Ehas, J. H., & Huchra, J. P. 1979, PASP, 91, 35. 
Whitmore, B., Gilmore, D., & Jones, C. 1993, ApJ, 407, 489. 

Zwicky, E., Herzog, E., Wild, P., Karpowicz, M., & Kowal, C. T. 1968, Catalogue of Galaxies 
and Clusters of Galaxies^ California Institute of Technology, Pasadena. 



This manuscript was prepared with the AAS E^TjrjX macros v3.0. 



- 39 - 



Figure Captions 

Fig.l — A schematic diagram of the configuration of the 4-Shooter focal plane during the 
gathering of the survey data. The four CCDs are denoted as 1, 2, 3, and 4. The arrow on 
the left is the CCD readout direction (along the columns) and is the motion of the sky across 
the detectors. 

Fig. 2 — The response curves of the 4-Shooter F555W (V4) and F785LP [I4) filters used in 
the survey. The response includes the telescope, CCD, 4-Shooter camera, and filter quantum 
efficiencies. The upper curve for each filter is the response with no atmosphere; the lower 
curve is the response at an airmass of f .2 at Palomar. 

Fig. 3 — (a) The k-correction as a function of redshift and morphology for the V4 filter, (b) 
The k-correction as a function of redshift and morphology for the I4 filter, (c) The V4—I4 
color as a function of redshift and morphology. 

Fig. 4 — The photometric offset as a function of time for four pairs of adjacent scans. The 
upper two panels show offsets for scans taken on the same night; the bottom two panels 
show offsets for scans taken more than a month apart. The time is derived from the object 
declination, which is linearly related to the time elapsed since the start of the scan. 

Fig. 5 — (a) The differential galaxy surface density as a function of V4 isophotal magnitude 
for the 5 survey fields, (b) The same for the I4 band. 

Fig. 6 — Residuals between the mean galaxy surface density for all five fields and the galaxy 
surface density in each individual field as functions of V4 and I4 isophotal magnitude. Solid 
line represents mean residual level, dashed line gives the prediction based on the observed 
angular two point correlation function (Neushaefer et al. f99f). 

Fig. 7 — (a) The galaxy-to-star count ratio as a function of V4 isophotal magnitude for 
all objects in the f3h field (solid dots) and for those objects which lie within 200 < 
CCD Column < 600. The results for the other survey fields are similar, (b) The same 
for the I 4 images. 

Fig. 8 — (a) V4 band star counts as a function of CCD column number. The I4 band results 
are similar, (b) V4 band galaxy counts as a function of CCD column number. The I4 band 
results are similar. 

Fig. 9 — (a) The V4 fiux filter as a function of redshift shown in intervals of O.f over the 
z range [0.2,f.2]. The ordinate is the amplitude of the filter. Numbers indicate the filter 
redshift. No color evolution is assumed, (b) The I4 fiux filter. 

Fig.fO — The radial filter as a function of redshift. The ordinate is the amplitude of the 
filter. A core radius of fOO/i~^ kpc and cutoff radius of lh~^ Mpc are used. The 5 radial 
filters shown, from broadest to narrowest, are for z = 0.2, 0.3, 0.4, 0.6, and f .0. 

Fig.ff — The histogram of S{i^j) values generated using a matched filter tuned to detect 
clusters at z = 0.6 in the 0^ field. 
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Fig. 12 — The cluster signal, S'c/, as a function of characteristic magnitude m* and passband. 
The corresponding filter redshift is shown along the upper abscissa. Points for a cluster at a 
given redshift are connect by solid lines. The solid point for each cluster denotes the redshift 
estimate that the matched filter algorithm would assign to the cluster. The estimates shown 
here are derived prior to the application of the cluster signal correction. The vertical dashed 
lines indicate the survey magnitude limits. 

Fig. 13 — The corrected cluster signal, S^i^ ^ as a function of characteristic magnitude m* 
and passband. The corresponding filter redshift is shown along the upper abscissa. Points 
for a cluster at a given redshift are connect by solid lines. The solid point for each cluster 
denotes the redshift estimate that the matched filter algorithm would assign to the cluster. 
Note the improved redshift estimate accuracy provided by the application of the CSC. 

Fig. 14 — The offset between estimated redshift and true redshift, Zest — -^true, as a function of 
true redshift for ~ 400 simulated clusters. Solid points represent clusters which are detected 
at the > 4(7 level. Crosses represent clusters detected below the 4(7 level. 

Fig. 15 — (a) The offset between estimated redshift and true redshift, Zest — Ztme^ as a function 
of the corrected V4 band matched filter signal, S^f^ = {S{i^j) — Sbg)/C SC , for 6 real clusters, 
(b) The same for the /4 band. 

Fig. 16 — The match-filtered /4 band galaxy catalog as a function of filter redshift for the 
region centered on the GHO cluster 1322+3028 {zohs = 0.70; cluster 063 in Table 4). The 
filter redshift is listed at the bottom of each 14.5' X 14.5' subimage; north is at the top and 
west to the right. The GHO cluster 1322+3027 {zobs = 0.75; cluster 059 in Table 4) is the 
peak to the east of the central cluster. 

Fig. 17 — The Abell richness Nr[< 1.0h~^ Mpc) as a function of for the PDCS clusters 
and for simulated clusters in the V4 and I4 passbands. The solid dots, open squares, and 
stars indicate cluster redshift estimates in the range 0.2 < Zest ^ 0.4, 0.5 < Zest ^ 0.7, and 
Zest ^0.8, respectively. 

Fig. 18 — Left: the mean spurious detection surface density in 10 simulated galaxy catalogs 
as a function of size and peak significance. Right: the mean surface density of observed 
cluster candidates as a function of size and peak significance. 

Fig. 19 — The matched filter richness, Ac/, as a function of redshift, richness, k-correction, 
and passband for the 3520 simulated clusters with r"^'"* profiles. The dashed lines show 
the relationship for clusters with richness class 1-4 and with elliptical galaxy k-corrections. 
The solid lines show the relationship for clusters with Sbc k-corrections. The points are the 
actual detected clusters (open circles are supplemental clusters). The shaded area represents 
region below the typical 3(7 detection threshold. Some matched cluster candidates lie within 
the shaded area in one passband but the significance of the detection in the other passband 
is above the 3(7 limit. 

Fig. 20 — The selection function for clusters with r"^'"* profiles. The dashed lines show the 
function for clusters with richness class 1-4 and with elliptical galaxy k-corrections. The 
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solid lines show the lunction lor clusters with Sbc k-corrections. II 4 separate solid lines are 
not seen it is because the selection probabilities for richer clusters are identical. 

Fig. 21 — The cumulative number of clusters per square deg function of redshift and 

richness for the PDCS, the GHO catalog, the Abell/ACO catalog, the APM cluster catalog, 
and the Edinburgh-Durham-Milano (EDM) cluster redshift catalog. The upper plot shows 
PDCS cluster counts as a function of the Abell richness estimator. Nr. The lower plot shows 
PDCS cluster counts as a function of the richness estimator. For clarity, the APM and 
FDM values are shown on separate plots. 

Fig. 22 — The differential redshift distributions for clusters detected in the V4 and I4 bands. 
The curves show the expected distributions for clusters consisting of non-evolving elliptical 
(solid line) and Sbc (dashed line) galaxies for Oo = 0.2 and Oo = 1- The Oo = 1 curves 
always predict fewer clusters at a given redshift. See text for further details. 

Fig. 23 — The cumulative redshift distributions for clusters detected in the V4 and I4 bands. 
The curves show the expected distributions for clusters consisting of non-evolving elliptical 
(solid line) and Sbc (dashed line) galaxies for Oo = 0.2 and Oo = 1- The Oo = 1 curves 
always predict fewer clusters at a given redshift. See text for further details. 
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Plate Captions 

Plate 1 — Finding charts for PDCS clusters 001 — 020. The images are of the central 
134" X 134" regions; north is up and west is to the right. The cluster ID is given at the 
top of each picture. The clusters are shown in the passband in which the significance of the 
peak signal is highest. Occasionally there are narrow gaps in the images - these are just the 
unusable between 4-Shooter CCDs. 

Plate 2 — Finding charts for PDCS clusters 021 - 040. 

Plate 3 — Finding charts for PDCS clusters 041 — 060. 

Plate 4 — Finding charts for PDCS clusters 061 — 079 and supplemental PDCS cluster 001 
(supplemental clusters have an S after the ID number). 

Plate 5 — Finding charts for supplemental PDCS clusters 002 — 021. 

Plate 6 — Finding charts for supplemental PDCS clusters 022 — 028. 



TABLE 1. PDCS Fields 



Center Coordinates EfF. Area (deg^) 

R.A. (J2000) Dec. / b Ay, A^ h 



OQh 29m IP +05° 31' 55" 112.76 -56.9 0.055 0.027 1.074 0.970 

02^^ 28"^ 33^ +00° 55' 45" 166.64 -53.5 0.028 0.014 1.120 0.927 

09^ 52"^ 57^ +47° 33' 20" 170.21 +49.8 0.004 0.002 1.052 0.995 

13^ 26"^ OP +29° 52' 55" 54.89 +81.9 0.031 0.015 1.062 1.014 

16^ 06"^ 17^ +41° 32' 40" 65.92 +47.9 0.019 0.009 1.088 1.026 
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TABLE 2. PDCS Photometric Standards 



Standard V4 I4 



HD 19445 8.04 7.77 

Ross 34 11.10 9.70 

Feige 67 11.71 12.70 

Ross 484 10.81 9.93 

BD +26 2606 9.71 9.46 

BD +28 4211 10.40 11.39 

BD +17 4708 9.45 9.20 
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TABLE 3. Matched Filter and Cluster Detection Parameters 

Parameter Value 

M^^ [h = 0.75, qo = 0.5) -21.00 

Ml (h = 0.75, qo = 0.5) -21.90 

a (LF slope) -1.10 

Tc (Core Radius) 100h~^ kpc 

r-co (Cutoff Radius) IM'^ Mpc 

Detection Threshold 95*^ pctl 

Peak Signal (S''^"*^) (dual band detection) max(Sy^"'^, Sj^^"*^) > Sa^g 

Peak Signal (single band detection) > Aabg 

V4 magnitude limit 23.8 

I4 magnitude limit 22.5 



TABLE 4. The Palomar Distant Cluster Catalog 



R.A. Dec. Va h 
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TABLE 5. The Supplemental Palomar Distant Cluster Catalog 



R.A. Dec. Va h 



ID 


(J2000) 




(J2000) 


Rad 




Peak 


(J 






N„ 






Rad 




Peak 


(J 


Ac; 




N„ 






001 


00 


30 


18 


.9 


+05 


46 


30 








.50 





.26 


1 


.86 


34 


.7 


22+ 
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63+ 9 
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.5 
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.9 


+05 


55 


05 








.40 





.39 


2 


.21 


31 
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NOTES TO TABLES 4 AND 5 
Cluster 005 : GHO cluster 0025+0455. 

Clusters 009 and 013 : These two detections may be part of the same system. Cluster 013 
not detected separate system in the I4 band. 

Cluster Oil : GHO cluster 0027+0535 (also known as 00 aa). 

Cluster 012 : GHO cluster 0025+0539. 

Cluster 017 : GHO cluster 0227+0024 (also known as 02 KPa). 

Clusters 026 and 027 : These clusters are not detected in the I4 band. Low V4 richness values. 

Cluster 028 : I4 band data for this cluster is not available. It would have been detected in 
the I4 band otherwise. 

Cluster 029 : Cluster was not detected in the I4 band. Low V4 richness values. 
Clusters 030 and 045 : These two detections may be associated with one another. 
Cluster 032 : GHO cluster 0949+4732. 

Cluster 033 : This detection may be associated with cluster 032. 
Cluster 038 : GHO cluster 0947+4758. 
Cluster 041 : GHO cluster 0951+4813. 

Clusters 043 : Cluster was not detected in the I4 band. Low V4 richness values. 

Clusters 044 : Cluster was not detected in the I4 band. Low V4 richness values. Eaint galaxies 
in V4 data lie below I4 completeness limit. 

Cluster 045 : Cluster was not detected in the V4 band. Near bright star. 
Cluster 049 : Abell 1739. 

Clusters 050 : GHO cluster 1324+2948. This detection may be associated with Abell 1739 
(Cluster 049). 

Cluster 051 : GHO cluster 1321+2949. 

Cluster 052 : GHO cluster 1323+2953 (also known as 13 KPN/3). 
Cluster 055 : GHO cluster 1324+3014 (also known as 13 KPA). 
Cluster 056 : GHO cluster 1324+3015. 
Cluster 058 : GHO cluster 1325+3021. 

Cluster 059 : GHO cluster 1322+3027 (also known as 13 KPNa). 
Cluster 063 : GHO cluster 1322+3028 (also known as 13 KPN7). 

Cluster 064 : Cluster was not detected in the V4 band. Eaint galaxies near 4-Shooter pyramid 
shadow region. 

Cluster 073 : GHO cluster 1604+4144 (also known as 16 7). 

Cluster 078 : Cluster was not detected in I4 band. Eaint galaxies in V4 data lie below I4 
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completeness limit. High V4 richness values. 
Cluster 079 : Cluster was not detected in V4 band. 

Supplemental Cluster 017 : This detection may be associated with Abell 1739 (cluster 049). 
Supplemental Cluster 020 : GHO cluster 1322+3014. 

Supplemental Cluster 023 : Background cluster behind GHO cluster 1605+4119. 
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Fig. 1. — A schematic diagram of the configuration of the 4-Shooter focal plane during the 
gathering of the survey data. The four CCDs are denoted as 1, 2, 3, and 4. The arrow on 
the left is the CCD readout direction (along the columns) and is the motion of the sky across 
the detectors. 
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Fig. 2. — The response curves of the 4-Shooter F555W (V4) and F785LP [I4) filters used in 
the survey. The response includes the telescope, CCD, 4-Shooter camera, and filter quantum 
efficiencies. The upper curve for each filter is the response with no atmosphere; the lower 
curve is the response at an airmass of 1.2 at Palomar. 




Fig. 3. — (a) The k-correction as a function of redshift and morphology for the V4 filter. 
The k-correction as a function of redshift and morphology for the I4 filter, (c) The V4 
color as a function of redshift and morphology. 
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Fig. 4. — The photometric offset as a function of time for four pairs of adjacent scans. The 
upper two panels show offsets for scans taken on the same night; the bottom two panels 
show offsets for scans taken more than a month apart. The time is derived from the object 
declination, which is linearly related to the time elapsed since the start of the scan. 
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. 5. — (a) The differential galaxy surface density as a function of V4 isopliotal magnitude 
the 5 survey fields, (b) The same for the I4 band. 
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Fig. 6. — Residuals between the mean galaxy surface density for all five fields and the galaxy 
surface density in each individual field as functions of V4 and I4 isophotal magnitude. Solid 
line represents mean residual level, dashed line gives the prediction based on the observed 
angular two point correlation function (Neushaefer et al. f99f). 
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Fig. 7. — (a) The galaxy-to-star count ratio as a function of isophotal magnitude for 
all objects in the 13h field (solid dots) and for those objects which lie within 200 < 
CCD Column < 600. The results for the other survey fields are similar, (b) The same 
for the /4 images. 



13hr Field: V band data 
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Fig. 8. — (a) V4 band star counts as a function of CCD column number. The I4 band results 
are similar, (b) V4 band galaxy counts as a function of CCD column number. The I4 band 
results are similar. 
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Fig. 9. — (a) The V4 flux fllter as a function of redshift shown in intervals of 0.1 over the 
z range [0.2,1.2]. The ordinate is the amphtude of the fllter. Numbers indicate the fllter 
redshift. No color evolution is assumed, (b) The I4 flux fllter. 
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Fig. 10. — The radial filter as a function of redshift. The ordinate is the amplitude of the 
filter. A core radius of 100h~^ kpc and cutoff radius of lh~^ Mpc are used. The 5 radial 
filters shown, from broadest to narrowest, are for z = 0.2, 0.3, 0.4, 0.6, and 1.0. 




Fig. 11. — The histogram of S{i^j) values generated using a matched filter tuned to detect 
clusters at 2; = 0.6 in the 0^ field. 




Fig. 12. — The cluster signal, S'c/, as a function of characteristic magnitude m* and passband. 
The corresponding filter redshift is shown along the upper abscissa. Points for a cluster at a 
given redshift are connect by solid lines. The solid point for each cluster denotes the redshift 
estimate that the matched filter algorithm would assign to the cluster. The estimates shown 
here are derived prior to the application of the cluster signal correction. The vertical dashed 
lines indicate the survey magnitude limits. 




Fig. 13. — The corrected cluster signal, S^i^^ as a function of characteristic magnitude m* 
and passband. The corresponding filter redshift is shown along the upper abscissa. Points 
for a cluster at a given redshift are connect by solid lines. The solid point for each cluster 
denotes the redshift estimate that the matched filter algorithm would assign to the cluster. 
Note the improved redshift estimate accuracy provided by the application of the CSC. 
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Fig. 14. — The offset between estimated redshift and true redshift, Zest — Ztme^ as a function of 
true redsliift for ~ 400 simulated clusters. Solid points represent clusters which are detected 
at the > 4(7 level. Crosses represent clusters detected below the 4(7 level. 




Fig. 15. — (a) The offset between estimated redshift and true redshift, Zest — -^true, as a 
function of tlie corrected V4 band matclied filter signal, S^i^ = {S{i^j) — Sbg) /C SC ^ for 6 
real clusters, (b) The same for the I4 band. 




Fig. 16. — The match-filtered I4 band galaxy catalog as a function of filter redshift for the 
region centered on the GHO cluster 1322+3028 {zobs = 0.70; cluster 063 in Table 4). The 
filter redshift is listed at the bottom of each 14.5' X 14.5' subimage; north is at the top and 
west to the right. The GHO cluster 1322+3027 {zobs = 0.75; cluster 059 in Table 4) is the 
peak to the east of the central cluster. 
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Fig. 17. — The Abell richness Nr[< 1.0h~^ Mpc) as a function of for the PDCS clusters 
and for simulated clusters in the V4 and I4 passbands. The solid dots, open squares, and 
stars indicate cluster redshift estimates in the range 0.2 < Zest ^ 0.4, 0.5 < Zest ^ 0.7, and 
Zest ^0.8, respectively. 
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Fig. 18. — Left: the mean spurious detection surface density in 10 simulated galaxy catalogs 
as a function of size and peak significance. Right: the mean surface density of observed 
cluster candidates as a function of size and peak significance. 




Fig. 19. — The matched filter richness, Ac/, as a lunction ol redshilt, richness, k-correction, 
and passband lor the 3520 simulated clusters with r"^'"* profiles. The dashed lines show 
the relationship for clusters with richness class 1-4 and with elliptical galaxy k-corrections. 
The solid lines show the relationship for clusters with Sbc k-corrections. The points are the 
actual detected clusters (open circles are supplemental clusters). The shaded area represents 
region below the typical 3(7 detection threshold. Some matched cluster candidates lie within 
the shaded area in one passband but the significance of the detection in the other passband 
is above the 3(7 limit. 




Fig. 20. — The selection function for clusters with r"^'"* profiles. The dashed lines show the 
function for clusters with richness class 1-4 and with elliptical galaxy k-corrections. The 
solid lines show the function for clusters with Sbc k-corrections. If 4 separate solid lines are 
not seen it is because the selection probabilities for richer clusters are identical. 
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Fig. 21. — The cumulative number ol clusters per square deg function of redshift and 

richness for the PDCS, the GHO catalog, the Abell/ACO catalog, the APM cluster catalog, 
and the Edinburgh-Durham-Milano (EDM) cluster redshift catalog. The upper plot shows 
PDCS cluster counts as a function of the Abell richness estimator. Nr. The lower plot shows 
PDCS cluster counts as a function of the richness estimator. For clarity, the APM and 
FDM values are shown on separate plots. 




Fig. 22. — The differential redshift distributions for clusters detected in the V4 and I4 bands. 
The curves show the expected distributions for clusters consisting of non-evolving elliptical 
(solid line) and Sbc (dashed line) galaxies for Oo = 0.2 and Oo = 1- The Oo = 1 curves 
always predict fewer clusters at a given redshift. See text for further details. 
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Fig. 23. — The cumulative redshift distributions for clusters detected in the V4 and I4 bands. 
The curves show the expected distributions for clusters consisting of non-evolving elliptical 
(solid line) and Sbc (dashed line) galaxies for Oo = 0.2 and Oo = 1- The Oo = 1 curves 
always predict fewer clusters at a given redshift. See text for further details. 



